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Application  of  Solidification  Theory  to 
Rapid  Solidification  Processing 


1.  Technical  Report  Summary 

This  semi-annual  technical  report  for  ARPA  Order  3751  covers  the 
period  October  1,  1979  to  March  31,  1980.  Work  is  reported  in  the  following 
areas  of  rapid  solidification  processing:  (1)  Interface  Stability  during 
Rapid  Solidification,  (2)  Eutectic  Solidification  and  the  Formation  of 
Metallic  Glasses,  and  (3)  Thermodynamics  of  Metastable  Equilibria. 

Manuscripts  have  been  prepared  for  publication  in  each  of  these  areas  and 
will  be  published  in  the  Proceedings  of  the  Second  International  Conference 
on  Rapid  Solidification  Processing:  Principles  and  Technologies.  This 
conference  was  held  March  24-26,  1980  in  Reston,  Va.  These  papers  are  included 
here  as  the  main  body  of  this  semi-annual  report.  In  addition,  a  paper 
entitled  Rapid  Solidification,  which  reviews  earlier  work  here,  is  Included 
as  a  final  item  in  this  report.  This  paper  was  presented  at  the  conference 
on  Laser  and  Electron  Beam  Processing  of  Materials  held  November  1979  in 
Boston,  Massachusetts,  and  will  be  published  in  the  proceedings  of  that 
meeting. 

Task  Objective 

The  objective  of  this  work  is  to  investigate  the  theory  of  rapid  solidi¬ 
fication  to  determine  whether  major  aspects  of  rapid  solidification  processes 
can  be  explained  by  conventional  solidification  theory  and  to  examine  non- 
equilibrium  effects  which  can  arise.  In  particular,  segregation  effects, 
glass-forming  tendencies,  and  rules  governing  the  formation  of  equilibrium 
and  non-equilibrium  phases,  including  alloy  composition  limits,  reaction 
sequences  and  metastable  phase  formation,  will  be  investigated. 


Technical  Problem  and  General  Methodoloi 

During  the  freezing  of  alloys,  a  planar  solid- liquid  interface  can 
become  unstable;  this  leads  to  a  cellular  or  dendritic  interface  and  causes 
solute  microsegregation.  The  transition  between  planar  and  non-planar 
growth  during  the  directional  solidification  of  binary  alloys  is  being 
investigated  using  morphological  stability  theory.  Conventional  constitutional 
supercooling  theory  predicts  that  the  stability  of  planar  solidification 
interfaces  will  decrease  as  solidification  velocities  increase.  At  large 
solidification  velocities,  however,  morphological  stability  considerations 
predict  that  planar  interface  stability  will  increase  with  solidification 
velocity.  Calculations  are  being  done  for  rapid  solidification  conditions  to 
determine  the  effect  of  solidification  velocity  on  interface  stability. 

Eutectic  alloy  compositions  are  frequently  found  to  provide  favorable 
conditions  for  formation  of  metallic  glasses.  The  theory  of  how  rapidly  one 
needs  to  solidify  to  obtain  a  glassy  structure  has  usually  focussed  on  the 
nucleation  of  the  solid  phase  as  the  difficult  step  in  crystallization.  In 
the  present  work,  the  focus  instead  is  on  the  growth  process.  Dlffusional 
sorting  of  alloy  components  is  necessary  to  form  a  two-phase  eutectic  com¬ 
posite  structure.  Theoretical  predictions  can  be  made  of  the  limitation  this 
imposes  on  eutectic  spacing  and  at  sufficiently  rapid  solidification  rates 
this  dlffusional  sorting  would  no  longer  be  possible,  perhaps  resulting  in 
amorphous  solids  being  formed.  Experimental  methods  designed  to  solidify 
eutectic  alloys  at  known  velocities  and  temperature  gradients  and  to  test 
this  theory  by  providing  quantitative  correlation  of  alloy  microstructure, 
such  as  eutectic  spacing  or  glass  formation,  with  solidification  conditions 
are  being  developed.  Two  methods  are  being  pursued:  (1)  rapid  directional 
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solidification  at  controlled  speeds  up  to  10  cm/s  and  (2)  electron  beam 
surface  melting  for  solidification  at  higher  velocities. 


Rapid  solidification  processes  are  frequently  regarded  as  leading  to 
non-equilibrium  conditions  which  cannot  be  predicted  by  theories  restricted 
to  conventional  equilibrium  thermodynamics .  Nevertheless,  a  careful  con¬ 
sideration  of  thermodynamic  phase  reactions  and  limits  on  metastable  phases 
may  adequately  describe  a  number  of  rapid  solidification  effects.  The 
thermodynamic  theory  of  multiphase  systems  is  being  investigated  to  determine 
rules  that  should  control  metastable  reaction  sequences,  relative  stability 
of  heterogeneous  microstructure  and  compositional  limits  imposed  by  the 
thermodynamics  during  rapid  solidification. 

Technical  Results,  Important  Findings  and  Conclusions 

Detailed  methods  and  results  are  described  in  the  main  body  of  the 
report.  Important  new  results,  findings  and  conclusions  include: 

-  (a)  The  phenomenon  of  solid-liquid  Interface  stability  during  directional 

solidification  of  a  binary  alloy  was  reexamined  with  special  emphasis  on  very 
rapid  solidification  rates.  Although  at  low  to  moderate  solidification  rates 
an  increase  in  solidification  velocity  leads  to  a  decrease  in  interface 
stability,  it  was  found  at  very  rapid  solidification  rates  (above  0.01  m/s  in 
copper  containing  aluminum) ,  the  perturbation  theory  predicts  a  strong  increase 
in  stabilization.  There  appear  to  be  two  dominant  effects  leading  to  this 
result.  The  first  effect,  known  as  absolute  stability  can  occur  even  if 
there  is  local  equilibrium  at  the  solid-liquid  Interface.  This  effect  arises 
because  only  short  wavelength  perturbations  are  important  at  high  velocities 
and  these  are  stabilized  by  surface  energy.  Numerical  calculations  for 
aluminum  with  various  concentrations  of  copper  and  a  liquid  temperature 


gradient  of  2(10  )  K/m  show  that  absolute  stability  is  important  for  solidification 
velocities  above  about  0.01  m/s;  in  fact,  for  copper  concentrations  less  than 

_3 

6x10  wt.  %,  the  interface  is  stable  at  any  solidification  rate.  The  second 
effect  is  caused  by  departure  from  local  equilibrium  at  the  solid-liquid 
interface.  Although  the  exact  forms  of  such  departures,  let  alone  their 
magnitude,  is  not  well  known,  the  perturbation  theory  can  be  modified  to 
include  their  salient  features.  Most  departures  from  local  equilibrium  lead 
to  increased  stability.  For  example,  all  constitutional  effects  will  vanish 
as  the  distribution  coefficient  approaches  unity. 

(b)  The  relationship  between  eutectic  solidification  and  the  ease  of 
formation  of  metallic  glasses  is  being  investigated.  For  many  systems, 
crystallization,  including  partitionless  crystallization,  of  alloys  into  a 
single  phase  solid  is  impossible  over  a  wide  range  of  composition  near  stable 
or  metastable  eutectics.  This  fact  forces  alloys  to  crystallize  into  two- 
phase  solids.  Because  of  the  need  for  diffusional  sorting  of  the  components 
apd  creation  of  solid-solid  surfaces,  the  kinetics  of  eutectic  crystallization 
are  relatively  slow  and  may  be  closely  related  to  the  ease  of  glass  formation 
of  these  alloys. 

Experiments  were  done  onthe  directional  solidification  of  Pd-6  at  % 

Cu-17  at  %  Si  alloys  which  show  the  evolution  of  microstructure  as  a  function 
of  interface  velocity.  A  structure  of  dendrites  and  interdendritic  eutectic 
at  low  velocity  0.25  ram/s)  becomes  a  fine  eutectic-like  structure  at 
intermediate  velocities  (y  1  mm/s)  and  finally  the  alloy  forms  glass  when 
crystallization  is  attempted  at  velocities  greater  than  2.5  mm/s. 

(c)  The  thermodynamics  of  metastable  phase  equilibrium  can  be  shown  to 
be  closely  related  to  that  of  stable  phase  equilibrium.  Metastable  phase 


equilibrium  often  can  be  represented  by  extensions  of  the  curves  on  phase 
diagrams  that  represent  stable  phase  equilibria.  Thermodynamic  principles 
rank  phase  equilibria  but  not  phases  in  a  hierarchy  of  increasing  stability. 
Thermodynamics  imposes  precursor  rules  in  multicomponent  systems  in  which  a 
prior  reaction  is  required  before  a  phase  can  appear.  A  given  hierarchy 
persists  over  a  domain  on  the  phase  diagram  bounded  by  surfaces  on  which  the 
ranking  of  two  equilibria  changes.  Multicomponent  phases  per  se  do  not  form 
a  hierarchy  because  a  given  phase  can  disappear  and  reappear  in  a  sequence  of 
spontaneous  processes.  Processing  for  producing  metastable  phases  must  place 
the  system  within  a  domain  in  which  the  desired  phases  can  form  spontaneously 
from  the  available  phases. 

(d)  Modifications  were  made  to  an  electron  beam  melting  apparatus  at 
NBS  to  allow  improved  surface  melting  experiments  to  be  performed. 

Special  Comments 

NBS  personnel  from  this  contract  played  major  roles  in  organizing  the 
Second  International  Conference  on  Rapid  Solidification  Processing  held  March 
24-26,  1980  in  Reston,  Va.  Attendance  at  the  conference  was  more  than  170. 
The  Chairman  of  the  Conference  Steering  Committee  was  Robert  Mehrabian,  Chief 
of  the  NBS  Metallurgy  Division  and  the  Conference  Secretary-Treasurer  was 
John  Manning.  This  conference  was  co-sponsored  by  DARPA,  ONR,  AFML,  ARO  and 
NBS.  A  conference  proceedings,  containing  more  than  40  papers,  will  be 
published.  Robert  Mehrabian  is  Chairman  of  the  Conference  Publications 
Committee.  Three  of  the  papers  written  for  this  proceedings  by  NBS  authors, 
which  report  work  done  on  the  current  contract,  are  included  as  part  of  this 
semi-annual  report. 


During  the  period  covered  by  this  semi-annual  report.  Dr.  R.  Schaefer, 


formerly  of  the  Naval  Research  Laboratory,  joined  our  staff  and  has  begun 

work  on  electron  beam  melting  aspects  of  this  work.  Prof.  T.  Massalski  qf _ 

Carnegie-Mellon  University  also  is  at. NBS  in  the  Metallurgy  Division  this _ 

.year  jsmd^  is  contributing"  as  a  consultant  to  the  current  work. 

Plans  and  Implications  for  Future  Research 

Since  the  current  work  has  shown  that  glass  formation  can  be  obtained 
during  directional  solidification  of  Cu-Si-Pd  alloys,  this  system  will  be 
studied  in  more  detail  to  permit  a  semi-quantitative  evaluation  of  equations 
concerning  the  effect  of  solidification  velocities.  Future  work  will  include 
attempts  to  measure  interface  temperature,  temperature  gradient,  and  the 
eutectic  spacing  along  with  an  identification  of  phases  present  during 
eutectic  solidification  at  velocities  just  below  the  transition  to  the  glass. 
Additionally,  experiments  on  the  electron  beam  melting  of  Al-Ag  and  Al-Cu 
alloys  will  be  initiated  to  test  morphological  stability  theory  at  high 

s' 

velocity  and  to  determine  the  critical  conditions  for  the  onset  of  massive 
(partitionless)  solidification,  respectively.  .  .  . 

Future  work  on  theoretical  aspects  of  this  problem  will  include 
investigation  of  solute  segregation  due  to  curved  solid-liquid  interfaces, 
the  effect  of  non-constant  velocities  on  Interface  stability  during  rapid 
solidification,  and  further  calculations  of  interface  stability  during 
growth  into  supercooled  liquids.  Equations  will  be  developed  to  bring  into 
morphological  stability  theory  the  crystal  anisotropies  and  apply  the  results 
to  high  solidification  velocities.  The  tendency  that  cells  may  have  to  follow 
the  heat  flow  direction  rather  than  to  take  on  crystallographic  forms 
characteristic  of  ordinary  cells  and  dendrites  will  be  Investigated. 
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2.  Report  of  Technical  Progress  and  Results 
Four  manuscripts  were  prepared  for  publication  during  the  period  covered 
by  this  semi-annual  report.  Since  these  papers  present  the  major  results 
obtained  in  work  on  this  contract  they  are  included  here  as  the  report  on 
technical  progress  and  results.  The  papers  on  "Interface  Stability  During 
Rapid  Solidification,"  "Eutectic  Solidification  and  Formation  of  Metallic 
Glasses,"  and  "Thermodynamics  of  Metastable  Equilibria"  will  be  published  in 
the  Proceedings  of  the  Second  International  Conference  on  Rapid  Solidification 
Processing,  held  March  1980  in  Reston,  Virginia.  The  paper  on  "  Rapid 
Solidification"  will  be  published  in  Laser  and  Electron  Beam  Processing  of 
Materials,  edited  by  C.  W.  White  and  P.  S.  Peercy,  Materials  Research  Society 
(proceedings  of  conference  held  Nov.  1979  in  Boston,  Massachusetts). 


Interface  Stability  During  Rapid  Solidification 


S.  R.  Coriell 

national  Bureau  of  Standards 
Washington,  DC  20234 

and 

R.  F.  Sekerka* 

Carnegie-Kellon  University 
Pittsburgh,  PA  15213 


ABSTRACT 

The  phenomenon  of  solid-liquid  interface  stability  during  directional 
solidification  of  a  binary  alloy  is  reexamined  with  special  emphasis  on  very 
rapid  solidification  rates.  For  ordinary  solidification  rates,  the  predictions 
of  the  perturbation  theory  of  morphological  stability  lead  to  results  that  are 
similar  to  those  implied  by  constitutional  supercooling;  however,  at  very  rapid 
solidification  rates,  the  perturbation  theory  predicts  a  vast  increase  in 
stabilization  in  comparison  to  constitutional  supercooling.  There  appear  to  be 
two  dominant  effects.  The  first  effect  can  occur  even  if  there  is  local 
equilibrium  at  the  solid-liquid  interface;  it  is  known  as  absolute  stability 
and  arises  because  only  short  wavelength  perturbations  are  important  at  high 
velocities  and  these  are  stabilized  by  surface  energy.  Numerical  calculations 
for  aluminum  with  various  concentrations  of  copper  and  a  liquid  temperature 
gradient  of  IflO1*)  K/m  show  that  absolute  stability  is  important  for  solidifi¬ 
cation  velocities  above  about  O.Cu  m/s;  in  fact,  for  copper  concentrations  less 
than  6xl0-3  wt.  *,  the  interface  is  stable  at  any  solidif ication  rate.  The 
second  effect  is  caused  by  depai-ture  from  local  equilibrium  at  the  solid-liquid 
interface.  Although  the  exact  forms  of  such  departures,  let  alone  their 
magnitude,  is  not  well  known,  the  perturbation  theory  can  be  modified  to  in¬ 
clude  their  salient  features.  Most  departures  from  local  equilibrium  lead  to 
increased  stability.  For  example,  all  constitutional  effects  will  vanish  as 
the  distribution  coefficient  approaches  unity.  Finally,  other  factors  are 
examined  with  an  aim  toward  the  identification  of  a  rea2m  where  experimental 
test  of  the  theory  would  be  meaningful. 


*  Consultant,  National  Bureau  of  Standards 
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I.  Introduction 


For  about  the  last  thirty  years,  the  principle  of  constitutional  super- 

CD 

cooling  has  been  used  as  a  guideline  to  ascertain  the  growth  conditions  that 

result  in  solid-liquid  interface  shape  instability  during  alloy  solidification. 

For  about  the  last  twenty  of  those  thirty  years,  this  principle  has  been 

(2-4) 

complemented  and  extended  by  the  theory  of  morphological  stability.  The 

approaches  of  these  analyses  are  somewhat  different,  constitutional  super¬ 
cooling  being  largely  an  application  of  thermodynamics  to  decide  if  a  liquid  is 
supercooled  (and  therefore  deemed  to  be  unstable)  with  due  respect  to  its 
composition,  and  morphological  stability  theory  being  a  detailed  kinetic  analysis 

of  the  growth  or  decay  of  a  perturbation  according  to  the  laws  of  heat  flow  and 

(5-9) 

diffusion.  Extensive  reviews  are  available  m  the  literature. 

Despite  the  different  approaches  of  the  constitutional  supercooling  principle 
and  morphological  stability  theory,  there  are  a  large  number  of  common  situations 
where  they  lead  to  similar  results,  or  at  least  similar  trends  within  the 
inaccuracies  of  our  knowledge  of  system  and  material  parameters.  For  this 
reason,  the  constitutional  supercooling  principle  is  widely  used  because  of  its 
simplicity.  Nevertheless,  there  are  situations  v?here  these  two  approaches  can 
lead  to  quite  different  results.  One  such  situation  occurs  when  the  average 
temperature  gradient  (namely  (k^g+k^G^)/ (ks+k^) ,  where  Gg  and  G^  are  temper¬ 
ature  gradients  in  solid  and  liquid,  respectively,  and  kg  and  k^  are  the 
corresponding  thermal  conductivities)  is  significantly  different  from  G^,  the 
only  temperature  gradient  that  the  constitutional  supercooling  principle 
considers.  A  second  more  dramatic  situation  occurs  in  the  case  of  melting,  as 
analysed  by  -  :.n  and  Jackson, for  which  constitutional  superheating  under¬ 

estimates  stability  by  orders  of  magnitude  because  it  uses  the  diffusivity  of 
the  solid  rather  than  the  diffusivity  of  the  liquid  (where  the  relevant  solute 
transport  occurs) . 

A  third  situation  where  constitutional  supercooling  and  morphological 
stability  lead  to  quite  different  answers  -  indeed,- the  main  subject  of  the 
present  paper  -  is  the  case  of  rapid  solidification.  At  very  rapid  rates  of 
solidification,  the  perturbation  theory  predicts  a  vast  increase  in  stabiliza¬ 
tion  in  comparison  to  constitutional  supercooling.  There  appear  to  be  two 
dominant  effects.  The  first  effect  can  occur  ever,  if  there  is  local  equilibrium 
at  the  solid-liquid  interface;  it  is  known  as  absolute  stability^  and  arises 
because  only  short  wavelength  perturbations  are  important  at  high  velocities 
and  these  are  stabilized  by  surface  energy.  This  case  will  be  presented  in 
Section  II  and  illustrated  with  specific  application  to  dilute  alloys  of  Cu  in 
Al.  The  second  effect  is  caused  by  departure  from  local  equilibrium  at  the 


solid-liquid  interface  and  is  taken  up  in  Section  III.  Section  IV  is  devoted 
to  the  special  case  of  a  supercooled  liquid  and  the  results  applied  to  the 
solidification  of  small  droplets. 

We  emphasize  that  the  analyses  of  the  present  paper  deal  with  constant 
solidification  rates  and  initially  planar  interfaces,  the  case  of  non-constant 
high  solidification  rates  -  although  typical  of  many  experimental  situations 
-  has  only  been  tractable  in  special  cases. 

II.  Local  Equilibrium  and  the  Transition  to  Absolute  Stability 

We  consider  the  morphological  stability  of  a  planar  solid-liquid  interface 

during  unidirectional  solidification  of  a  binary  alloy  at  constant  velocity  V. 
(12) 

We  outline  a  standard  linear  time  dependent  stability  analysis  and  then 
focus  our  attention  to  large  velocities  where  absolute  stability  prevails  and 
our  results  differ  significantly  from  those  of  constitutional  supercooling. 

We  choose  an  (x,y,z)  coordinate  system  (moving  with  the  planar  interface) 
such  that  the  solid-liquid  interface  is  described  by  z  =  W(x,y,t)  where  c  is 
the  time.  We  assume  that  W(x,y,t)  and  any  of  its  derivatives  are  sufficiently 
small  that  any  nonlinear  terms  can  be  neglected. 

We  solve  the  differential  equations 

(3TL/3t)  -  +  V(3Tl/3z),  (la) 

(3T  /3t)  -  k  V2T  +  V(3T  /3z) ,  (lb) 

5  S  S  S 

,  j 

(3c/3t)  =  DVzc  +  V(3c/3z) ,  (lc) 

where  T.  and  T  are  temperatures  in  the  liquid  and  solid,  respectively,  c  is 

L  S 

the  concentration  of  solute  in  the  liquid  (diffusion  in  the  solid  is  neglected), 
and  <s  are  thermal  dif fusivities  of  liquid  and  solid,  respectively,  and  D  is 
the  diffusion  coefficient  of  solute  in  the  liquid. 

The  boundary  conditions  far  from  the  solid-liquid  interface  are 


(3Tl/3z) 

»  exp 

(-Vz/c^) 

(2a) 

(3Ts/3z) 

-  Gs  exp 

<-Vz/,-8) 

(2b) 

c 

-  C 

•m 

z  *  «, 

(2c) 

-10- 


vhcrc  the  temperature  gradients  and  Ca 
constants.  Wc  consider  G  >  0  and  G,  >  0; 
Section  IV. 


and  the  solute  concentration  arc 
the  case  <  0  will  be  discussed  in 


The  boundary  conditions  at  the  solid-liquid  interface  are 


tl-t. 

(3a) 

V  -  V  +  (3W/3t) 

(3b) 

V  -  (k./Lv)(3Ts/3z)  - 

(kj/L^OTj/Dz) 

(3c) 

v  ■  -  D(3c/3z)/(c-c  ) 
s 

(3d) 

T  -  T  2  T.  -  T  FK  + 
h  e  ii  M 

g(c) 

(3c) 

cs  “  Ve)* 

(3f) 

where  v  is  the  interface  velocity,  and  k^  are  thermal  conductivities  of 
solid  and  liquid,  respectively,  L v  is  the  latent  heat  per  unit  volume,  cg  is 
the  concentration  in  the  solid  at  the  interface,  Tf[  is  the  melting  point  of  a 
flat  interface  in  the  absence  of  solute,  T  is  the  ratio  of  the  solid-liquid 
surface  tension  y  to  the  latent  heat  per  unit  volume,  K  «*  -  (3zW/3xz)  - 
(3zW/3y2)  is  the  interface  curvature,  the  dependence  of  melting  point  on 
solute  is  given  by  the  function  g(c)  and  the  relationship  between  c#  and  c  is 
given  by  the  function  ho(c).  Equations  (3e)  and  (3f)  arc  based  on  the 
assumption  of  local  equilibrium  at  the  solid-liquid  interface.  In  Equation 
(3c)  we  have  ommitted  correct  ions  that  might  arise  from  consideration  of 
surface  entropy  and  also  a  possible  dependence  of  g(c)  and  hg(c)  on  curvature. 


In  carrying  out  a  linear  stability  analysis,  wc  write  the  temperature  and 

concentration  fields  as  a  sum  of  an  unperturbed  part,  which  is  a  function  of  z 

alone,  and  a  perturbed  part  which  is  of  the  form  F  (z)  exp  [ot+i(m  x+u  y) ] , 

j  3c  y 

where  is  a  function  of  z  alone.  The  perturbed  solid-liquid  interface  is 
given  by 


*  ■  V(x,y,t)  -  4  exp  [ot  +  i(u  x  +  w  y)J,  (4) 

*  7 

where  4  is  the  perturbation  amplitude  at  t  ■  0,  and  and  w ,  are  spatial 

frequencies.  The  interface  is  unstable  if  the  real  pare  of  the  time  constant 

o  i s  positive  for  any  perturbation,  i.c.,  any  real  values  of  w  and  u  .  The 

x  y 

interface  is  stable  If  the  real  part  of  j  is  negative  for  all  perturbations. 


(3) 
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Upon  solving  the  differencial  equations,  one  finds 

C  -  V{- Wl“V/'£L>  “  W0s+VAs>  -  2kTMr“2“  + 

2taaC  ci(a-V/D)/  (a-pV/D)  }/(L  V  +  2kmC  a/(o-pV/D)}  (5) 

C  VC 

wlch 

a  -  (V/2D)  +  [(V/2D)2  +  w2  +  o/D],^ 

-  (V/2kl)  +  [(V/2kl)2  +  u2  +  o/kJ*5 

o  -  -(V/2.C  )  +  [  (V/2<  )2  +  u2  +  o/k  ]** 

S.'S  s  s 

•  *  (Vs  +  VL)/(2k)’ 

where  2k  -  k  +  k_  ,  u2  -  oi2  +  to2  G  -  Vc.(k-1)/D,  k  *  c  / cn,  m  ■  (dg/dc) 
sis  x  y  c  u  *  u 

and  p  -  1  -  (dh^/dc)  with  the  derivatives  evaluated  at  c  *  and 
defined  by  the  equation  *  hg(Cg) ,  i.e. ,  is  the  solute  concentration  in 
the  liquid  at  the  planar  solid- liquid  interface.  If  one  assumes  that  the 
functions  g(c)  and  hQ(c)  are  linear  in  c,  then  m,  k,  and  p  “  1  -  k  are  con¬ 
stants  independent  of  c. 

For  a  given  alloy  and  processing  conditions,  we  can  calculate  o(w)  from 
equation  (5).  If  che  real  part  of  a  is  positive  for  any  value  of  a,  the 
interface  is  unstable.  Alternately  we  can  specify  the  real  part  of  o,  of,  and 

calculate  some  other  variable  such  as  c  as  a  function  of  u.  For  o  "0,  the 

•  r 

minimum  value  of  c  (u) ,  which  occurs  at  some  specific  value  of  u,  gives  the 

stability- instability  demarcation.  Except  where  otherwise  specified,  we  have 

also  assumed  that  the  imaginary  part  of  a,  a^,  is  zero  at  the  stability- 

instability  demarcation,  consistent  with  a  detailed  analysis  in  the  thermal 

(121 

steady  state  approximation. 

We  have  carried  out  numerical  calculations  for  alloys  of  aluminum  containing 
copper  for  the  following  values  of  material  parameters: 

k^  -  90. 7J/ (msK) ,  k8  -  210J/(msK), 

Ly  -  1.0vSCJ/m2,  T}ir  -  1.0(10"7)mK, 

-  4.2(10"S)m?/s,  r.a  -  8. l<10"5)ur/s, 

D  -  5<l0*°)m'/*.  k  -  0.14,  n  -  -2.5  K/(v\.  2). 
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The  curve  in  Fig.  X,  based  on  equation  (5),  shows  the  critical  copper  concen¬ 
tration  (below  which  the  planar  interface  is  stable)  as  a  function  of  growth 
velocity  V  for  a  temperature  gradient  in  the  liquid,  G^,  of  2(10**)  K/m.  The 
dashed  straight  line  with  negative  slope  is  the  constitutional  supercooling 
criterion,  mC  «  C,  ,  and  the  straight  line  with  positive  slope  is  the  absolute 
stability  criterion,  mG  -  kT..r(V/D)2(soe  eq.  (7)  et  soq.).  The  constitutional 
supercooling  criterion  is  a  good  approximation  at  low  velocities  and  the 
absolute  stability  criterion  is  a  good  approximation  at  high  velocities.  In 
Fig.  2,  the  wavelength  \  -  2n/u  that  corresponds  to  the  onset  of  instability,  l.e. 
to  the  value  of  id  for  which  c ^(w) .  is  a  minimum,  is  shown  to  be  a  monotonically 
decreasing  function  of  interface  velocity.  Although  results  are  given  in 
Figs.  1 -and  2  for  V  »  10  m/s,  the  theory  should  not  be  taken  seriously  at 
these  velocities  since  the  diffusion  length  D/V  is  of  atomic  dimensions. 

From  Fig.  2,  one  finds  that  (V/2<  )  <<  u  and  (V/2kt)  <<  w,  which  permits 

S  u 

the  use  of  the  thermal  steady  state  approximation  according  to  which  we  let 
the  thermal  dif fusivities  approach  infinity.  For  the  stability-instability 
demarcation  (a  *  0) ,  Equation  (5)  Chen  reduces  to 

0  -  Vu{-kLGL-k3Gs-2kTMro.2  +  2kmGc  (a-V/D)  /  (a-pV/D)  }/ 

(LyV  +  2kmCc<n/  (a-pV/D)  >.  (6) 

The  denominator  is  always  positive  so  that  only  the  properties  of  the  numerator 

(14) 

need  be  analyzed.  The  stability  criterion  may  be  written  as 

2^  +  VLv  >  2fanCcS(A,k),  (7) 

where  S  is  a  function  of  the  dimensionless  variables  k  and  A  •  kT,,rv2/(D2cG  ). 

tl  c 

For  A  _>  1,  the  interface  is  always  stable,  while  for  A  <  1,  the  spatial 
frequency  at  the  onset  of  instability  and  the  function  S(A,k),  which  lies 
between  zero  and  unity,  can  be  calculated  by  finding  the  one  and  only  real 
positive  root  of  the  cubic  equation 

r3  +  (2k- 1)  r  -  ^k/A*5)  -  0.  (8) 

The  quantities  S  and  u  are  then  given  by 

a 

S(A,k)  -  1  +  (A/4k)  (1  -  r2  +  2k  ’)  -  (3A>5r/2),  (9) 

u  -  (V/2D)  (r1*  -  l)*5. 

a 


(10) 


Rewriting  equation  (10)  in  the  form  r"  -  1  +  (2D10JV)2,  we  recognize  two 
limiting  eases,  viz.,  (2D<oJV)  »  1  and  (2VtoJV)  «  1.  First,  we  consider 
(2Dw  /V)  »  1,  and  consequently  r  >>  1.  The  solution  of  equation  (8)  is  then 
r3  =  2k/AJ*.  It  follows  that  S  *  1  -  (3k/r2)  and  to^  *  (V/2D)r2  a  [(kV/AkT^rD) 
(2kjCL  +  VLv)]173.  If,  in  addition,  VLv  «  21^0^,  then  o»B  *  [kk^G^/IkT^D]1^3, 
and  the  stability  criterion  is  (k^/k^  *  mGc,  which  is  shown  as  the  solid 
line  with  negative  slope  in  Fig.  1.  At  low  velocities  this  line  is  a  better 
approximation  to  the  curve  than  constitutional  supercooling  to  which  it  would 
be  identical  for  k  -  k_  .  The  proportionality^15^  between  to  and  V173  is  in 

S  L 

agreement  with  the  experimental  data  of  Morris  and  Winegard  on  lead- 
antimony  alloys. 

For  the  second  limiting  case,  (2Dw  /V)  <<  1  and  r  -*■  1.  If  we  let  r  *  1  +  e 

A 

with  0  <  c  «  1,  then 

c  =  [k/ (1+k)  ]  [  (l-A^J/A**]  , 

S  *  I(l-k/2)/  (1+k)  ]  (1-A*5)2, 

,  ( [kV2 (VL^kjC^  ] /  [  f  1+k)  (l-k/2) 2kTMrD2 ]  }% , 

and  the  stability  criterion  approaches  the  absolute  stability  criterion  A  “  1. 
For  VLv  »  2kG^ ,  we  have 

to  *  ( [kV3L  ]/[  (1+k)  (l-k/2) 2kTMrD2]}?l 

si  v  n 

Thus  for  small  velocities  (2Dub/V)  »  1  and  is  proportional  to  V173  while 
for  large  velocities  (2Doj  /V)  <<  1  and  is  proportional  to  v376.  These 
approximations  are  verified  by  the  results  in  Fig.  2.  Plots  of  S(A,k)  and 

(14) 

additional  analytic  results  are  available. 

Uc  close  this  section  with  two  additional  numerical  examples  for  Al-Cu 
alloys.  The  first  of  these  illustrates  the  effect  of  G^  on  the  critical 
concentration  for  instability  at  a  sufficiently  large  solidification  velocity 
that  constitutional  supercooling  is  a  poor  approximation.  From  Figs.  1  and  2, 
at  a  velocity  of  1  m/s  the  critical  concentration  is  0.21  wt.  Z  Cu  and  u  ■ 

t  m 

2.05  (107)m  Cor  G.'  ■  2(10ll)K/m.  Moreover  equation  (7)  shows  that  the 

L 

instability  criterion  iM.1 1  be  practically  independent  of  provided  that 
<<  VLv/2k^  ■  6(10f,)K/n-  Calculations  for  higher  values,  viz.  G^  »  106, 

107,  10?K/m,  give  critical  concentrations  of  0.22,  0.24,  and  0.34  wt.  «  and  u  ■ 

j  ® 

2.1(107),  2.6(107),  and  4.2(107)m”  ,  respectively,  indicating  a  relative 
insensitivity  to  the  liquid  temper ature  gradient. 
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The  second  numerical  example  illustrates  the  decree  to  which  an  instability 

can  develop  if  only  a  finite  time  is  available  before  solidification  is  complete. 

-4 

For  V  -  lm/s,  a  100  pm  layer  will  solidify  in  a  time  tfi  -  10  s.  In  order  to 
observe  an  unstable  interface,  the  instability  should  be  greatly  amplified  in 
a  fraction  of  this  time,  i.e. ,  we  require  Chat  oc  »  1  [see  eq.  (4)].  Taking 
t  -  O.lt  ■  10-^s,  we  require  o  >>  10s  s  Calculations  with  G.  *  2(104)K/m 
and  a  -  10*  and  107  s”  yield  critical  concentrations  of  0.23  and  0.36  wt.  X' 

and  ut  -  2.1(107)  and  3.4(107)b_1,  respectively.  These  concentrations  are  in 

n 

an  experimentally  accessible  range  but  the  wavelengths  involved  (M).2  um)  would 
necessitate  electron  microscopy. 

III.  Departure  from  Local  Equilibrium 

Ve  now  briefly  discuss  the  effect  of  departures  from  local  equilibrium  at 
the  solid-liquid  Interface  upon  interface  stability.  Ve  generalize  the  local 
equilibrium  boundary  conditions,  equations  (3e)  and  (3f),  to  the  form^ 


f(Te-Tx,  cr  Tx), 

(3e') 

-  h  (cx,  VTj), 

(3f') 

where  the  subscript  I  indicates  evaluation  at  the  solid-liquid  interface.  Ve 
require  that  v  -  0  and  h  -  h^  when  Te  ”  Tj  "  We  eacPand  the  above  equations 
In  a  Taylor  series  about  the  temperature  equilibrium  temperature  Tcq»  and 

solute  concentration  at  the  planar  solid-liquid  interface,  viz., 

v  -  V  +  PTl(Te-Tc0)  -  (Tj-Tjq) ]  +  Uc(cI-cI0>  +  WA(TrTI0), 

csl  ”  h^c10’Tc0'TI0)  +  kc(crcI0)  +  M(Te“’W  ~  (TI_TI0^ * 

where  uT  -  Sf/aCT^-Tj),  pc  -  3f/3cx,  uA  -  3f/3Tx>  kc  -  3h/3cx>  kj.  -  Sh/SC^-Tj) 
and  the  partial  derivatives  are  evaluated  at  T^q,  Txq,  and  c^^.  With  this 
generalization,  vc  obtain  in  place  of  equation  (5) 

o  -  VU-kjC^aj-V/^)  -ksG9(VV/<8)]UA  -  2^^  - 

+  2ta>'C  a(a-V/D)/(a-p*V/D))/[L  VU_ 

C  VI 

+  [2to'Cca/(a-p'V/D))UK},  (5') 

where  m*  »  a  +  (u  /u_),  p*  •  1  -  k  +  (k_ii  /u_) , 

Cl  C  1  C  T 
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UA  -  1-  [pa/pt3U  +  (kTa'V/D)/(c-p>V/D)]1 
UT  “  1  "  (yA/uT)  +  2k  a/avyTJ, 


and 


UR  -  1  +  lV2UT/<nTDGe)] ll  -  LvWA/(2k  a)]. 

As  p^  -*■  •,  we  approach  local  equilibrium  and  equation  (5')  approaches  equation 
(5). 

Although  equation  (5')  is  rather  general  in  allowing  for  deviations  from 
local  equilibrium,  our  knowledge  concerning  the  forms  of  the  functions  f  and  h 
is  extremely  limited  and  quantitative  conclusions  based  on  equation  (5')  are 
not  possible.  This  is  clearly  an  area  in  which  experimental  and  theoretical 
research  is  desirable.  We  will  attempt  to  draw  some  conclusions  in  a  few 
simplified  cases. 

Case  A  c  _  ■  c. 

-  si  I 

The  deviation  from  equilibrium  is  so  large  that  kc  “  1  and  Gc  "  0;  there 
are  no  constitutional  effects  and  the  situation  is  similar  to  solidification 
of  a  pure  material.  In  addition  to  there  being  no  instability,  there  is  no 
macrosegregation  in  contrast  to  the  stable  case  for  k£  f  i. 

Case  B  c  ,  »  h„(c)  and  v  ■  f(T  -T_) 

-  si  u  el 

Solid  composition  depends  only  on  liquid  composition  as  given  by  the  phase 
diagram  and  velocity  depends  only  on  deviation  of  the  temperature  T^  from  the 
temperature  T  .  We  have  p  »  p.  -  k_  “  0  so  that  m'  «  m,  p'  *  1  -  k  ,  U,  “ 

■  1,  and  »  1  +  2k  a/(Lvp^).  Thus,  the  numerator  of  equation  (S')  is  un¬ 
changed  from  the  local  equilibrium  result,  equation  (5).  Further,  the  denominator 
is  still  always  positive  since  p^  >  0;  the  stability- instability  demarcation 
is  unaffected  by  p^,.  However,  for  sufficiently  small  y  we  expect  that  a  will 
be  reduced  for  fixed  growth  conditions  or,  alternatively,  a  larger  value  of 
concentration  will  be  needed  to  obtain  a  given  value  of  o  at  fixed  V.  If 

PT  »  [2kaV]/(LvV  +  2kinCca/Ia-PV/D]}  =  D/(I '  r)  -  0.05m/sK 

holds  (where  the  simplification  to  D/T^r  is  valid  for  large  velocities),  the 
term  in  the  denominator  involving  Is  small  compared  to  the  sum  of  the  other 
two  terms,  and  has  little  effect.  Since  we  are  unaware  of  any  reliable 


-17- 


1 


measurements  of  y^  for  aluminum  alloys,  we  have  done  some  calculations  for 
values  of  y^  that  bracket  0.05  ra/sK.  These  calculations  based  on  equation 
(5')  are  summarized  in  Table  1;  it  follows  that  there  is  only  a  small  increase 
in  critical  concentrations  as  y^  decreases  from  •  to  0.1  m/sK  but  that 
critical  concentrations  nearly  double  as  y^  decreases  from  0.1  to  0.01  a/sK 
Note  that  if  f  is  a  linear  function  of  (T  -Tx),  then  a  velocity  of  1  m/s  and  a 
kinetic  coefficient  of  0.01  m/sK  correspond  to  a  kinetic  undercooling  of 
100K. 

Case  C  pure  material  and  v  »  f (T^-Tj,  Tj) 

Temkin  and  Polyakov^i9^  have  considered  a  kinetic  law  with  the  properties 
that  >  0  while  yT  -  yA  >  0  for  Tj  large  (near  Ty)  and  PT  -  VA  <  0  for  Tj 
sufficiently  small.  Such  a  kinetic  law  is  sketched  in  Fig.  3  for  two  different 
values  of  equilibrium  temperature,  T  ’  and  T^".  If  one  writes  v  »  F(Tft,  Tj) 
instead  of  v  *  f(T  -T.,  Tj) ,  then  yT  -  uA  -  -3F/3TX  and  yT  -  3F/3Tft,  permitting 
geometrical  interpretation  via  Fig.  3.  A  specific  analytic  example  is  v  “ 
U0(Te-Tx)  exp  (-O/Tj),  where  yQ  and  0  are  constants,  and  from  which  y^.  «  yQ 
exp  (-0/Tx)  and  bT  ~  UA  "  Ug  (-O/Tj)  [l-OC^-T^/T*] .  For  such  a  kinetic 
law,  equation  (5')  reduces  to 


a  -  ( ( VV/ Kj) -ksGs (as+V/<s)  ]  [ yr-uA]  -  y^kT^a}/ 

{Lv(yT-yA)  +  2k  a}.  (11) 

Note  that  y^  -  yA  ■  0  corresponds  to  the  maximum  in  velocity  as  a  function  of 

T_  for  fixed  T  . 

I  e 


So  long  as  y^  -  y^  >  0,  analysis  of  equation  (11)  leads  to  results  that  are 

qualitatively  similar  to  case  B.  If  y^  -  y^  <  0,  further  analysis  is  required 

since  the  sign  of  the  denominator  in  equation  (11)  is  not  obvious.  For 

—  W 

simplicity,  we  let  k  ■  k.  and  «  «  k  so  that  a  -  { (V/2ic  ) 2  +  ui2  +  o/«T  K 

S  L  S  L  L  L 


TABLE  I 


CRITICAL  CONCENTRATIONS  AND  SPATIAL  FREQUENCIES  FOR 
VARIOUS  VALUES  OF  0  AND  y_  AT  C.  -  Z(10")K/m  AND  V  -  1  m/s 


oCa"1) _ U^m/sK) 


0 


10s 

m 

10  7 

«B 

10s 

0.1 

107 

0.1 

106 

0.01 

107 

0.01 

„(wt.  2) 

O 

1 

Nl 

BE 

B 

0.21 

2. 05 

0.23 

2.1 

0.36 

3.4 

0.26 

2.9 

0.53 

5.5 

0.40 

4.S 

1.29 

9.4 
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Figure  3.  Sketch  of  interface  velocity  V  as  a  function  of  interface 
temperature  for  two  different  values  of  the  equilibrium  temperature. 

and 

o  =  {[2kgGga  -  VLvCa  -  V/2<L> 3 [uA-  VTJ  “  V^k^ru2®}/ 

(21^0  -  Lv(pa  -  uT)).  (12) 

We  consider  growth  into  a  supercooled  liquid  with  Gg  =  0.  It  is  easy  to  show 
Chat  a  *  0  is  not  a  solution  of  equation  (12)  for  any  value  of  u  when  (y^-y^) 

<  0.  This  suggests  that  we  need  to  allow  for  the  possibility  chat  o  is 
complex.  Numerical  solution  of  equation  (12)  indicates  that  this  is  the  case. 
Taking  <s  “  “  6.0(10-S)m2/s,  ■  kg  -  150J/(msK),  y^  »  O.Olm/sK,  y^  - 

U.  -  -0.01  m/sK,  L  »  1.08  CJ/m3,  T.J  -  1.0(10‘7)mK,  and  V  -  1  m/s,  we  find 

A  V  rl 

Chat  the  interface  is  unstable  for  <u  <  1. 6(104)  m"  ;  the  maximum  value  of  the 

real  part  of  a  occurs  at  u  =  8(103)  m  *  with  o  »  (2.9 (103 )  +  i  2.5(10'')]  s  1. 

Thus,  for  cither  sign  of  y^  -  y^,  the  interface  is  unstable  for  small  values 

of  u  and  is  stabilir.ed  by  surface  tension  for  large  values  of  u>;  however,  for 

(y„-y . )  <  0  the  nature  of  the  instability  is  oscillatory.  The  treatment  by 
1  A  ( 19 ) 

Temkln  and  Polyakov  did  not  uncover  these  oscillatory  instabilities  be¬ 
cause  it  used  steady  state  heat  flow  equations. 
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IV.  Ctowth  Into  Undercoolcd  Molts 


Solidification  into  undcrcooled  Belts  (GL<0)  deserves  special  attention 
because  for  this  case,  the  constitutional  supercooling  principle  always  predicts 
instability.  As  shown  below,  the  interface  can  be  stable  under  these  conditions 


for  sufficiently  large  values  of  Ca.  In  the  thermal  steady  state  approximation, 
the  previous  analysis,  equations  6  and  7,  is  still  valid  as  long  as  + 


k^  ■  2(4^  +  VLy  >  0;  however, 


if  this  inequality  does  not  hold,  equation 


(6)  predicts  instability  as  u  -*•  0.  This  suggests  that  a  more  rigorous  analysis 


is  required  since  the  thermal  steady  state  approximation  assumes  that  u  >> 


(V/2ica)  and  u  »  (V/2kl)  . 


We  let  o  -*■  0  and  define  the  quantity  R(w)  »  {-k^G^Ca^-V  /  k^)  -  ksGs(as+V/x:s) }/ 
{2lt^}.  If  R(w)  <  0,  the  thermal  field  is  stabilizing  while  for  R(w)  >  0  the 
thermal  field  is  destabilizing.  Clearly  R(0)  ”  -CgkgKL/(kk<s)  ant^  m 

(-kLGL-kgGs)/2k  and  these  have  opposite  signs  when  kjGL  +  kgGs  <  0.  These 
results  Indicate  that  the  thermal  field  is  stabilizing  for  very  small  y 
(provided  G^  >  0)  and  is  destabilizing  for  large  u  (when  k^C^  +  kaGg  <  0) . 


The  special  case  of  Gs  *  0  requires  a  more  delicate  analysis  which  we 
undertake  only  for  =•  kg  and  < ^  =  Ka>  Then  R(i»)  “  (-GL/2)(1  -  [1+(2icui/V)  Z]-*5} 
is  a  monotonically  increasing  function  of  w  with  R(0)  **  0  and  R(«)  **  (-G^/2)  >  0. 
Hence,  the  thermal  field  is  destabilizing  in  this  case.  Since  the  solute 
field  is  always  destabilizing  and  surface  tension  is  important  only  for  large 
u,  we  conclude  that  such  a  situation  is  unstable. 


In  applying  such  an  analysis  to  the  solidification  of  supercooled  droplets, 
there  is  an  additional  factor  that  needs  to  be  considered,  viz.,  that  in  a 
small  droplet  only  rather  large  values  of  oi  can  occur.  If  the  droplet  is  of 
radius  R,  then  the  largest  meaningful  perturbation  wavelength  1  "  2a/y  <_  2R 
ot  u  >  t/R.  Thus,  for  a  droplet  of  1  ym  diameter  w  _>  6(10e)m  1  so  that  w  > 
V/2<k  even  for  velocities  of  the  order  of  1  ra/s.  Thus,  R(w)  v  (-k^CL/2k) ,  and 
the  surface  tension  term  T,Xw2  >  3.6(10c)K/m.  For  G.  >  -1.2(107)K/m  or  V  < 

1.0  m/s,  we  would  expect  a  stable  interface  for  the  solidification  of  a 
droplet  of  pure  aluminum.  When  solute  is  present,  the  magnitude  of  the  surface 
tension  term  must  be  sufficiently  large  to  overcome  the  destabilizing  tempera¬ 
ture  and  solute  gradient  terms.  Numerical  calculations  with  G  «*  0  and  w  >_ 

r  -1  S  —4 

6(106)bi  yield  critical  concentrations  of  0.12,  0.02S,  and  7.7(10  )  wt.  X 
for  velocities  of  0.01,  0.1,  and  1,0  m/s,  respectively.  If  we  take  o/V  ■  106 
m  1  rather  than  zero,  we  find  critical  concentrations  of  0.15,  0.038,  and 
0.03'i  wt.  X  for  velocities  of  0.01,  0.1,  and  1.0  m/s,  respectively  with  y  * 
b(10r')  a  ns  United  by  the  particle  radius.  For  o/V  ■  107  m""*,  y  is 
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IP 


no  longer  Halted  by  the  particle  radius  and  we  find  critical  concentrations 
of  2.34,  0.37,  and  0.33  wt.  %  with  u  -  1.5(107),  1.7(107),  and  2.9(107)  sf1 

SI 

for  V  *>  0.01,  0.1,  and  1.0  m/s,  respectively,  Thus,  although  an  undercooled 
melt  with  Cg  »  0  is  unstable,  strictly  speaking,  for  any  composition,  some 
minimum  composition  is  needed  for  this  instability  to  be  observed  in  small 
droplets  because  of  capillary  stabilization. 

V.  Discussion 

We  have  demonstrated  in  some  detail  how  the  principle  of  constitutional 
supercooling  and  morphological  stability  theory  can  lead  to  quite  different 
results  at  rapid  solidification  rates.  In  particular,  this  occurs  because  of 
increased  capillary  stabilization  of  the  short-wavelength  perturbations  impor¬ 
tant  at  high  velocities  (absolute  stability) ,  because  of  increased  departures 
at  the  solid-liquid  interface  from  local  equilibrium,  and  because  of  the 
important  role  of  the  temperature  gradient  in  the  solid  in  the  case  of  under- 
cooled  melts.  We  have  also  presented  some  numerical  calculations  that  illustrate 
the  refinements  in  the  growth  conditions  needed  to  observe  instabilities  in 
the  cases  of  solidification  of  thin  films  and  fine  droplets.  These  effects 
suggest  the  need  for  more  careful  experimental  observations  at  rapid  solidifica¬ 
tion  rates. 

We  emphasize  that  all  of  the  above  conclusions  are  based  on  the  use  of 

macroscopic  transport  theory.  Surely,  when  solidification  rates  become  so 

rapid  that  critical  lengths,  such  as  D/V,  become  the  order  of  atomic  dimensions, 

the  use  of  macroscopic  transport  theory  cannot  be  justified.  For  Cu  in  Al,  D 
_q  .  • 

*  5  x  10  sr/s  and  D/V  would  be  less  than  1QA  for  V  >  5  m/s.  The  critical 
wavelength  for  instability  is  greater  than  D/V  for  large  V;  thus,  D/V  is  the 
length  with  which  we  have  to  be  concerned.  As  we  have  seen,  departures  from 
local  equilibrium  can  be  handled  in  a  formal  way  but  there  is  little  known 
about  exact  forms  and  magnitudes  of  these  departures.  Ve  caution  again  that 
our  results  are  only  strictly  applicable  to  constant  V;  the  case  of  large  buc 
time-dependent  V  is  probably  prevalent  in  most  experiments  that  involve  rapid 
solidification  rates  and  remains  as  a  ripe  subject  for  further  research. 
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ABSTRACT 

The  relationship  between  eutectic  solidification  and  the  ease  of  formation 
of  metallic  glasses  is  investigated.  For  many  systems,  crystallization, 
including  partitionless  crystallization,  of  alloys  into  a  single  phase  solid 
is  impossible  over  a  wide  range  of  composition  near  stable  or  metastable 
eutectics.  This  fact  forces  alloys  to  crystallize  into  two-phase  solids. 
Because  of  the  need  for  diffusional  sorting  of  the  components  and  creation  of 
solid-solid  surfaces,  the  kinetics  of  eutectic  crystallization  are  relatively 
slow  and  may  be  closely  related  to  the  ease  of  glass  formation  of  these  alloys. 

Experiments  are  reported  on  the  directional  solidification  of  Pd-6  at  Z  Cu- 
17  at  Z  Si  alloys  which  show  the  evolution  of  microstructure  as  a  function  of 
interface  velocity.  A  structure  of  dendrites  and  interdendritic  eutectic  at 
low  velocity  (<  0.25  mm/s)  becomes  a  fine  eutectic-like  structure  at  inter¬ 
mediate  velocities  1  mm/s)  and  finally  the  alloy  forms  glass  when  crystalli¬ 
zation  is  attempted  at  velocities  greater  than  2.5  mm/s. 
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Introduction 


Eutectic  systems  are  multicomponent  systems  in  which  liquid  remains  the 
equilibrium  phase  to  low  temperatures  because  no  solid  with  the  same  composition 
is  stable.  Because  of  the  stability  of  the  liquid  at  low  temperatures  these 
systems  have  been  prime  candidates  for  glass  formers  by  rapid  solidification. 

The  theory  and  practice  of  how  rapidly  one  needs  to  solidify  has  primarily 
focussed  on  the  nucleation  and  growth  of  a  single  solid  phase  as  the  difficult 
step  in  crystallization  (1-3).  We  have  focussed  instead  on  the  coupled 
growth  of  two  (or  more)  solid  phases. 

There  are  two  aspects  being  explored:  1)  the  thermodynamic  question  for 
each  system  of  the  absence  of  a  solid  phase  with  the  same  composition  as  the 
liquid  which  at  some  temperature  below  the  eutectic  becomes  stable  relative  to 
the  liquid,  2)  the  kinetic  question  of  how  fast  a  liquid  can  crystallize  into 
mixtures  of  crystals  of  differing  composition.  The  diffusional  sorting  of 
the  components  is  a  relatively  well  understood  slow  step  which  governs  the 
crystallization  rate.  These  ideas  have  led  to  a  series  of  predictions  for 
such  systems  which  are  being  explored  experimentally. 

The  thermodynamic  question  is  best  answered  by  examination  of  the  position 
of  the  Tq  curves  (A)  for  all  possible  solid  phases  on  the  phase  diagram.  A  T_ 
curve  is  the  intersection  of  the  solid  and  liquid  free  energy  curves  plotted 
on  a  phase  diagram.  It  marks  the  bound  on  compositions  of  a  solid  phase  which 
can  form  from  liquid  of  any  composition  at  that  temperature.  Figure  la  shows 
a  hypothetical  phase  diagram  with  a  stable  and  a  metastable  eutectic.  If  the 
composition  range  of  an  individual  phase  is  very  narrow  or  if  the  solidus  is 
retrograde  (5) ,  the  T^  curves  will  not  extend  far  from  that  phase  at  temperatures 
of  interest.  This  produces  a  rather  wide  range  of  alloy  compositions  between 
the  T_  curves  of  the  different  phases  which  must  crystallize  as  two-phase 
solids.  Even  partitionless  (massive)  crystallization  of  liquids  in  this 
composition  range  into  a  single  phase  is  impossible.  Wider  composition  ranges 
requiring  two  phase  crystallization  are  possible  if  a  particular  phase  fails 
to  nucleate  (for  example,  g  in  Figure  la) . 


(a) 

Figure  1.  a)  Hypothetic  phase  diagram  with  a  stable  and  a  metastable 
eutectic.  The  T_  curves  for  the  three  solid  phases  relative  to  the  liquid 
are  shown.  In  tne  shaded  region  of  composition  between  T«  curves, 
crystallization  to  a  single  phase  crystalline  solid  is  impossible, 
b)  Relationship  of  interface  temperature  to  interface  velocity  for 
coupled  growth  of  two  phases.  T^,  is  the  eutectic  temperature. 


The  growth  of  a  two-phase  solid  from  the  liquid  requires  the  diffusional 
sorting  of  the  components  in  the  liquid  phase  as  well  as  the  creation  of  new 
surface  area  between  solid  phases.  These  requirements  for  coupled  growth  lead 
to  well  known  relationships  (6,7)  between  interface  undercooling  AT  below  the 
eutectic  temperature,  interface  velocity  V  and  eutectic  spacing  X;  namely. 


and  redundantly 


D(T)A2 

“*T" 


(1) 

(2) 


A  =  2A2/AT.  (3) 

Here  D(T)  is  the  interdiffusion  coefficient  in  the  liquid  which  is  a  strong 
function  of  interface  temperature  at  large  undercoolings.  The  parameters  A^ 
and  A2  are  approximated  by  mAC/8  and  2y  /AS,  respectively,  where  n  is  an 
average  positive  liquidus  slope,  AC  is  the  composition  difference  between  the 
solid  phases,  y  is  the  a- 8  surface  energy  and  AS  is  an  average  entropy  of 
fusion  of  the  solid  phases.  The  above  equations  use  the  extremum  principle 
although  refinements  of  this  aspect  are  possible  (8).  Extension  of  eutectic 
solidification  theory  to  ternary  systems  has  been  performed  (9)  and  leads  to 
similar  results. 


For  a  given  interface  velocity,  the  interface  undercooling  required  for 
eutectic  solidification  based  on  equation  (1)  is  much  greater  than  for  single 
phase  crystallization.  In  fact  for  this  reason,  interface  attachment  kinetics 
are  neglected  in  equation  (1).  For  example,  for  velocities  of  10  pm/s  one 
expects  undercoolings  of  the  order  of  10  K  (10)  for  single  phase  crystallization 
of  metals,  whereas  for  eutectics  such  a  velocity  would  lead  to  an  undercooling 
of  the  order  of  IK  (11) .  For  larger  interface  velocities  the  temperature 
dependence  of  the  diffusion  coefficient  becomes  very  important  and  the  under¬ 
cooling  becomes  much  larger.  Figure  lb  shows  schematically  the  form  of  Eqn. 

1  for  a  diffusion  coefficient  with  Arrhenius  behavior.  A  maximum  velocity 
exists  beyond  which  eutectic  solidification  cannot  occur.  Such  a  maximum  also 
occurs  for  single  phase  solidification  but  it  occurs  at  a  much  higher  velocity. 
Similarly,  minimum  transformation  times  from  time-temperature  transformation 
curves  have  been  calculated  by  Boswell  &  Chadwick  (12) .  Alloys  required  to 
crystallize  at  velocities  greater  than  this  maximum  due  to  the  rate  of  heat 
extraction  and  which  cannot  crystallize  as  a  single  phase  material  will 
undoubtedly  form  a  glass. 


The  purpose  of  this  research  is  to  determine  whether  this  prediction  of  a 
maximum  crystallization  velocity  for  a  eutectic  is  related  to  the  ease  with 
which  certain  alloys  form  glass.  Preliminary  experiments  are  reported  in 
which  a  glass-forming  alloy  has  been  directionally  solidified  at  various 
interface  velocities.  Alloy  microstructure  including  the  transition  to  the 
glassy  state  is  described  as  a  function  of  interface  velocity. 


Experiments 

Rapid  directional  solidification  is  being  conducted  by  quenching  thin 
alloy-filled  quartz  tubes  (0.75  mm  I.D.  x  1.5  mm  O.D.  x  10  cm  long)  from  a 
furnace  at  1100  °C  into  liquid  Ga  at  various  controlled  speeds.  Samples  are 
initially  lowered  at  0.25  mm/s  to  obtain  about  3  cm  of  crystalline  solid. 

This  step  prevents  bulk  undercooling  of  the  liquid  and  permits  study  of  the 
growth  aspects  of  the  alloy  unencumbered  by  nucleation  difficulties.  Samples 
are  then  quenched  at  the  desired  speed.  An  important  consideration  here  is 
whether  the  interface  velocity  is  equal  to  the  quenching  rate.  Care  was  taken 
in  the  present  experiments  to  maximize  the  extraction  of  heat  from  the  samples 
during  quenching  by  the  use  of  thin  samples  with  thin  container  walls  and  a 
high  boiling  point  quenchant.  If  the  heat  flow  is  directional  and  parallel  to 
the  tube  axis  in  the  alloy,  the  Interface  velocity  must  be  less  than  or  equal 
to  the  quenching  rate.  In  fact  an  initial  transient  where  the  interface 
accelerates  to  the  quenching  speed  was  expected  and  observed  microstructurally 
to  be  about  2  cm  long.  -25- 


The  alloy  chosen  for  t ii6  research  Is  Pd-6  at  2  Cu-17  at  %  SI  because  of 
its  well  known  ease  of  glass  formation.  Unfortunately,  the  Pd-Si  as  well  as 
the  Pd-Cu-Si  phase  diagrams  are  not  well  known  (13-15) .  Figure  2a  and  2b  show 
micrographs  of  this  alloy  solidified  at  0.25  and  1.0  mm/s,  respectively.  Both 
samples  show  strong  alignment  of  microstructure  parallel  to  the  tube  axis  and 
indicate  directional  heat  flow  for  this  geometry  and  these  quenching  speeds. 

At  0.25  mm/s  (Fig.  2a)  the  microstructure  consists  of  both  faceted  and  unfaceted 
dendrites  and  an  interdendritic  two  phase  eutectic.  Suprisingly  both  types  of 
dendrites  appear  to  be  the  same  phase  and  have  a  composition,  determined  by 
electron  microprobe,  of  78  at  2  Pd,  20.6  at  %  Si  and  1.4  at  2  Cu  and  an 
electron  diffraction  pattern  which  has  been  indexed  based  on  an  orthorhombic 
cell  with  a  =  7.46X,  b  =  9.27X  and  c  =  9.60&.  The  phase  occurring  along  with 
this  phase  in  the  eutectic  has  a  composition  73  at  2  Pd,  12.1  a/o  Si  and  14.9 
a/o  Cu  and  appears  to  have  cubic  symmetry. 


Figure  2.  Longitudinal  sections 
of  Pd-Cu-Si  alloys  solidified  at 
various  velocities,  solidification 
direction  up;  (2)  0.25  mm/s, 
optical  micrograph;  (b)  1.0  mm/s, 
optical  micrograph;  (c)  transi¬ 
tion  from  dendritic  to  coupled 
growth  with  increasing  velocity, 
SEM  micrograph  (inverted  signal). 
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At  1.0  ran/s  the  dendritic  phase  is  absent  and  the  raicrostructure  consists 
of  a  fine  eutectic  structure  as  shown  in  Figure  2b.  A  complete  characterization 
of  this  structure  is  currently  underway.  The  transition  from  denuritic  to 
fine  eutectic  structure  with  increasing  velocity  is  quite  significant.  This 
transition  is  shown  dramatically  in  Figure  2c  which  occurred  in  the  initial 
transient  region  of  the  sample  quenched  at  1  irnn/s.  Similar  transitions  from 
dendritic  to  eutectic-like  structures  with  increasing  velocity  have  been 
observed  in  off-eutectic  Pb-Sn  alloys  by  Cline  and  Livingston  (16).  It 
appears  that  alloys  over  a  range  of  compositions  near  eutectics  can  undergo 
coupled  growth  of  two  or  more  phases  at  high  velocity.  Hence  the  kinetic 
treatment  presented  in  the  Introduction  applies. 

Alloys  quenched  at  2.5  mm/s  exhibit  a  very  different  microstructure  —  they 
are  amorphous.  In  the  initial  transient  there  is  the  transition  from  dendritic 
to  coupled  growth  already  described,  a  refinement  of  the  eutectic  spacing  and 
then  a  sharp  interface  as  shown  in  Figure  3a.  The  light  region  above  the 
sharp  interface  is  featureless  and  has  been  determined  to  be  glassy  from  its 
x-ray  and  electron  diffraction  patterns.  The  glassy  region  extends  the  re¬ 
maining  5  cm  length  of  the  sample  and  occupies  the  region  where  steady  state 
interface  speed  is  ordinarily  observed.  It  is  believed  that  the  general 
vicinity  of  the  sample  shown  in  Figure  3a  represents  the  place  where  the 
interface  velocity  reached  the  maximum  shown  in  Figure  lb.  The  transition 
from  the  dark  to  gray-etching  area  below  the  glass  interface  represents  a 
dramatic  reduction  in  the  eutectic  spacing.  However,  the  structure  appears  to 
remain  eutectic-like  right  up  to  the  glass  interface  as  seen  in  Figure  3b. 

This  region  of  fine  eutectic  spacing  may  possibly  be  the  material  that 
crystallized  on  the  part  of  the  kinetic  curve  (Fig.  la)  with  positive  slope. 

A  small  star-shaped  dendrite  visible  in  Figure  3a  formed  in  this  region, 
possibly  due  to  large  liquid  undercooling,  but  was  overtaken  by  the  eutectic 
interface. 


Figure  3.  Transition  from  fine  eutectic  growth  to  glassy  alloy  observed 
when  crystallization  at  2.5  mm/s  is  attempted,  optical  micrograph,  entire 
sample  width.  (b)  SEM  view  of  interface  between  fine  eutectic  structure 
and  glass.  -  -27- 
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A  conclusion  of  this  paper  is  that  the  maximum  velocity  of  eutectic 
solidification  of  this  alloy  is  between  1.0  and  2.5  mm/s.  Attempts  to  crystal¬ 
lize  at  velocities  greater  than  this  lead  to  glass  formation.  We  have  also 
seen  similar  sharp  interfaces  between  fine  eutectic  structures  and  glass  in 
electron-beam  surface  melted  samples. 

Future  work  will  include  attempts  to  measure  interface  temperature, 
temperature  gradient,  and  eutectic  spacing  along  with  an  identification  of  the 
phases  present  during  eutectic  solidification  at  velocities  just  below  the 
transition  to  the  glass.  This  information  will  permit  a  semi-quantitative 
evaluation  of  Eqn.  1. 
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ABSTRACT 

The  thermodynamics  of  metastable  phase  equilibrium  is  closely  related  to 
chat  of  stable  phase  equilibrium.  Metastable  phase  equilibrium  can  often  be 
represented  by  extensions  of  the  curves  on  phase  diagrams  that  represent 
stable  phase  equilibria.  Thermodynamic  principles  rank  phase  equilibria  but 
not  phases  in  a  hierarchy  of  increasing  stability.  Thermodynamics  imposes 
precursor  rules  in  multicomponent  systems  in  which  a  prior  reaction  is 
required  before  a  phase  can  appear.  A  given  hierarchy  persists  over  a  domain 
on  the  phase  diagram  bounded  by  surfaces  on  which  the  ranking  of  two  equilibria 
changes.  Multicomponent  phases  per  se  do  not  form  a  hierarchy  because  a  given 
phase  can  disappear  and  reappear  in  a  sequence  of  spontaneous  processes. 
Processing  for  producing  metastable  phases  must  place  the  system  within  a 
domain  in  which  the  desired  phases  can  form  spontaneously  from  the  available 
phases. 
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Introduction 


Thermodynamics  Is  the  science  of  the  Impossible.  It  enables  us  to  tell 
with  certainty  what  cannot  happen,  but  Is  noncommittal  about  the  things  that 
are  possible.  Thermodynamics  is  at  its  best  when  nothing  more  can  happen,  a 
condition  called  equilibrium.  The  concept  of  equilibrium  has  been  fruitfully 
extended  to  reversible  processes  which  are  on  the  verge  of  being  impossible. 

When  a  system  is  not  at  equilibrium,  thermodynamics  retains  its  certainty 
about  the  impossible.  However,  the  domain  of  the  thermodynamically  possible 
may  be  so  large  that  thermodynamics  has  lost  much  of  its  usefulness  as  an  aid 
to  prediction.  In  a  complicated  nonequilibrium  system  such  as  we  often 
encounter  in  rapid  solidification,  thermodynamics  would  be  almost  totally 
useless  if  applied  to  the  entire  system,  for  then  an  irreversible  process  in 
one  part,  combined  with  an  impossible  process  in  another,  could  lead  to  an 
overall  reduction  in  the  appropriate  thermodynamic  potential  and  a  verdict  of 
"possible".  Hence,  it  is  a  common  assumption  that  we  apply  thermodynamics 
locally,  unless  we  know  specifically  that  there  are  long-range  interactions. 
Similarly,  we  often  assume  that  different  processes  are  Independent,  and  apply 
thermodynamic  criteria  to  each  individually. 

Another  important  assumption  is  that  of  constrained  equilibrium,  in  which 
certain  processes  are  assumed  not  to  occur  or  to  occur  so  slowly  that  thermo¬ 
dynamics  is  fully  applicable  to  all  the  other  proccssses  remaining.  These 
three  assumptions  permit  us  to  apply  thermodynamics  locally  and  to  individual 
processes,  but  the  validity  of  such  assumptions  must  be  checked  experimentally. 

An  example  of  constrained  equilibrium  is  metastable  phase  equilibrium. 

Here,  the  constraint  is  that  one  or  more  of  the  stable  phases  is  absent.  The 
phases  that  are  actually  present  reach  equilibrium  subject  to  this  constraint. 
This  definition  of  metastable  phase  equilibrium  focusses  on  the  absence  of  one 
or  more  of  the  equilibrium  phases,  and  does  not  ascribe  any  unusual  charac¬ 
teristic  to  the  metastable  ones  present. 

A  second  example  of  constrained  phase  equilibrium  occurs  when  the  solidifi¬ 
cation  is  so  rapid  that  one  or  more  of  the  components  cannot  redistribute 
among  the  phases  in  the  time  scale  of  the  experiment.  In  a  completely 
partitionless  solidification  to  a  single  solid  phase,  the  temperature  where 
the  liquid  and  solid  have  equal  free  energies  traces  out  a  surface  (T  )  on 
the  multicomponent  phase  diagram  between  the  liquldus  and  solidus  for°that 
phase  and  its  extrapolation.  Below  this  surface,  provided  diffusion  is 
limited,  it  becomes  thermodynamically  possible  to  solidify  a  liquid  completely 
even  though  the  system  may  be  above  its  solidus. 

In  this  paper  we  will  focus  on  metasfable  equilibria  and  partitionless 
(sometines  called  massive)  transformations,  their  representation  on  phase 
diagrams  and  the  hierarchical  laws  governing  the  thermodynamically  possible 
sequences  of  phases  and  raetastable  phase  equilibria.  Certain  basic  concepts 
and  relationships  described  in  earlier  reviews  (1,2)  and  standard  textbooks 
will  be  assumed,.  We  begin  with  the  single  component  case  and  highlight  the 
major  complications  introduced  for  multicomponent  systems.” 

Metastable  Phase  Equilibria  in  Single  Component  Systems 

A  single  phase  can  reach  a  mecastable  phase  equilibrium  when  nothing  more 
is  possible  except  the  formation  of  new  phases.  Nucleation  theory  tells  us 
that  there  are  natural  barriers  to  the  formation  of  new  phases  from  mecastable 
ones.  The  range  of  conditions  under  which  a  phase  can  be  metastable  is 
bounded  by  kinetic  factors.  Imposed  by  nucleation  and  growth,  and  thermodynamic 
limits  cocsnonly  called  spinodals. 
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What  is  important  is  that  it  is  reasonable  to  assume  that  there  are  no 
discontinuities  in  thermodynamic  properties  as  a  phase  moves  from  stable  to 
metastable.  The  same  thermodynamic  measurements  that  are  made  on  stable 
phases  can  be,  and  have  been,  made  on  metastable  phases.  Heat  capacity, 
volume,  vapor  pressure,  and  other  such  properties  are  rigorously  defined  and 
have  been  measured.  From  such  measurements,  energy,  entropy,  and  free  energy 
can  be  rigorously  determined.  Metastable  phases  obey  the  usual  solution  laws. 
If  chcy  are  dilute,  Henry's  and  Raoult’s  laws  apply.  They  can  even  be  ideal: 
to  wit,  a  supersaturated  vapor  at  a  density  far  below  the  critical-point 
density. 

In  a  single  component  system  phases  are  in  equilibrium  when  the  chemical 
potentials  p  of  the  component ,  which  in  this  case  are  the  same  as  the  molar 
Gibbs  free  energies,  in  the  several  phases  are  equal.  A  graph  of  the  free 
energies  of  several  phases  as  a  function  of  temperature  at  constant  pressure 
is  shown  in  Figure  1.  At  each  temperature  the  stable  phases  are  the  ones  with 
the  lowest  value  of  F.  Consistent  with  the  phase  rule  only  one  phase  is 
stable  except  where  the  two  lowest  curves  cross.  All  other  phases  are 
motastable.  With  changing  pressure  the  free  energy  curves  will  shift  at  a 
rate  proportional  to  the  volumes  of  each  phase  and  different  phases  may  become 
stable. 
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Figure  1.  Chemical  potential  (free  energies)  of  various  one- 
component  stable  and  metastable  phases  as  a  function  of 
temperature  at  constant  pressure.  The  stable  phase  at  each 
temperature  is  the  one  with  the  lowest  p.  Stable  and  metastable 
two-phase  equilibria  occur  where  two  curves  cross.  Between 
crossings  the  hierarchy  of  phases  ranked  according  to  their 
values  of  p  is  listed. 

A  graph  of  the  domains  where  each  phase  is  stable  is  called  a  phase  diagram 
and  is  shown  in  Figure  2.  It  consists  of  areas  in  which  a  phase  is  stable 
bounded  by  curves  in  which  two  phases  are  the  equilibrium.  Three  curves  come 
together  at  triple  points  of  three-phase  equilibrium.  Vapor-liquid  critical 
points  where  a  two-phase  curve  terminates  are  of  little  concern  in  rapid 
solidification.  The  axes  can  be  any  two  independent  combinations  of  T,  P,  or  p 
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Figure  2.  A  one— component  phase  diagram  with  T.  P.  or  u  as  axes*  showing 
stable  two  and  three-phase  equilibria  as  curves  and  (triple)  points. 

Equilibria  between  metastable  phases  can  also  be  mapped  on  phase  diagrams. 

Each  of  the  curves  in  the  metastable  phase  diagram  represents  a  crossing  of 
two  free  energies.  As  in  the  case  of  the  stable  diagram  each  curve  is  identified 
with  the  equilibrium  of  the  two  phases,  one  of  which  is  stable  on  one  side  of 
the  curve,  the  other  is  stable  on  the  other.  In  Figure  3  this  sidedness  of 
each  curve  has  been  identified  by  labelling  the  stabler  phase  on  each  side. 

The  properties  of  these  two-phase  equilibria  are  unaffected  by  whether  or  not 
there  exists  a  more  stable  phase.  If  in  some  part  of  the  diagram  the  two- 
phase  equilibrium  is  stable,  the  curve  representing  it  can  be  extrapolated 
beyond  the  triple  point  where  it  becomes  metastable.  The  stable  triple  points 
also  have  become  the  crossings  of  three  two-phase  curves.  The  six  rays 
emanating  from  a  stable  triple  point  are  alternately  stable  and  metastable. 


Figure  3.  A  metastable  phase  diagram  with  the  same  phases  as  shown  in  Figure 
2.  Triple  points  have  become  triple  crossings  with  metastable  equilibria 
indicated  by  light  lines.  Within  each  region  a  given  phase  hierarchy  persists. 
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Such  a  metastable  phase  diagram  also  features  points  where  three  metastable 
phases  are  in  equilibrium  and  where  three  two-phase  curves  cross. 

The  metastable  diagram  is  seen  to  consist  of  many  regions  bounded  by  no 
more  than  three  curves.  Within  each  region  there  is  a  definite  hierarchy  of 
phases  according  to  their  free  energies  which  is  unchanged  over  the  entire 
region.  At  a  bounding  curve  a  pair  of  phases  change  place  in  the  hierarchy 
and  a  different  hierarchy  exists  in  each  adjacent  region.  In  region  X  the 
phase  hierarchy  i»  a  ■»  t  +  8  +  while  in  region  Y  the  order  of  6  and  y  have 
been  reversed  by  the  traversing  of  the  metastable  B  -  y  equilibrium  curve. 

Each  triple  point  represents  a  place  there  three  phases  have  equal  p. 

Three  two-phase  equilibria  curves  representing  the  three  pairings  of  the  three 
phases  must  cross  there.  The  crossings  of  two  curves  representing  two  pairs 
of  two  phases  with  none  in  common  are  not  triple  points. 

The  thermodynamic  hierarchy  for  a  one-component  system  ranks  phases  by 
their  molar  free  energies.  In  a  thermodynamic  processing  path  of  heating, 
cooling,  and  pressure  changes,  a  phase  can  be  retained  and  moved  into  a  region 
where  it  is  metastable  and  above  one  or  more  phases  in  the  hierarchy.  During 
spontaneous  phase  changes  it  can  only  move  down  the  hierarchy.  A  liquid 
undercooled  just  below  its  equilibrium  melting  curve  becomes  metastable  only 
with  respect  to  the  stable  crystal.  It  does  not  become  metastable  with 
respect  to  a  metastable  crystal  until  it  has  been  brought  below  the  melting 
curve  of  the  metastable  phase.  This  metastable  melting  curve  can  often  be 
observed  as  a  reversible  phenomenon,  and  can  often  be  estimated  from  the 
stable  phase  diagram  if  the  metastable  phase  becomes  stable  elsewhere  in  phase 
diagrams. 

Figure  3  included  only  phases  which  were  stable  in  some  portion  of  Figure 
2.  All  curves  in  Figure  3  are  extrapolations  of  curves  from  the  stable  phase 
diagram.  Metastable  equilibria  as  well  as  the  metastable  hierarchies  of  all 
the  phases  that  are  stable  in  some  part  of  a  stable  phase  diagram  are  thus 
easy  to  estimate.  For  example,  bismuth  can  be  made  to  crystallize  at  ambient 
pressure  into  a  high  pressure  phase  which  appears  to  have  a  metastable  melting 
point  of  175  *C  consistent  with  the  metastable  extensions  in  the  phase  diagram 
(3).  It  is  not  surprising  that  most  observed  metastable  phases  appear  as 
stable  phases  somewhere  in  phase  diagrams  either  at  different  pressure,  temper¬ 
ature  or  with  small  additions  of  other  components.  In  order  for  a  phase  to  be 
formed  by  rapid  solidification  it  has  to  be  within  a  few  hundreths  eV  of  the 
stable  phase  and  it  is  unlikely  that  a  phase  will  come  that  close  to  the 
lowest  free  energy  without  breaking  through  at  some  adjacent  pressure,  temperature 
or  alloy  addition.  It  is  not  impossible  though.  Figure  4  gives  a  hypothetical 
case  where  a  metastable  phase  6  comes  close  to  being  stable  without  ever 
becoming  so,  and  having  a  free  energy  curve  which  is  not  parallel  to  that  of 
any  other  phase.  It  is  quite  straightforward  to  construct  the  hypothetical 
free  energy  surfaces  corresponding  to  the  metastable  equilibria  of  Figure  4. 

While  there  is  no  singularity  in  thermodynamic  properties  of  the  liquid  at 
any  of  the  stable  or  metastable  melting  points,  there  is  a  rapid  change  in 
heat  capacity  near  the  glass  transition  (4).  The  heat  capacity  of  the  amorphous 
phase  at  the  glass  'transition  is  cooling  rate  and  heating  rate“dependent  but 
seems  at  the  slowest  rates  to  approach  a  well  behaved  reversible  limit  (5,6). 

For  the  faster  rates  there  is  a  lack  of  reversibility,  but  heating  and  cooling 
curves  can  establish  an  upper  and  lower  bound  to  the  free  energy  of  a  glass 
when  its  properties  depend  on  its  thermal  history.  The  slow  cooling  limit  of 
the  glass  transition  seems  well  defined,  and  indicates  that  the  free-energy 
curves  of  glass  and  liquid  merge  smoothly,  exhibiting  only  a  rapid  change  in 
curvature.  On  phase  diagrams,  the  curves  representing  the  metastable  equilibria 
between  the  amorphous  phase  and  some  other  phase  arc  continuous  through  the 
glass  transition,  but  exhibit  a  rapid  change  in  slope.  If  the  metastable 
crystalline  phase  does  not  catalyze  the  nucleation  of  a  stabler  crystalline 


phase,  the  reversible  growth  or  dissolution  of  this  phase  in  contact  with  a 
glass  should  be  observable. 
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Figure  4.  A  metastable  phase  diagram  similar  to  Figure  3  except  that 
the  5-phase  is  nowhere  stable. 


Solidification  processing  permits  a  climbing  in  free  energy  away  from  the 
equilibrium.  The  solidification  cycle  is  not  unlike  a  heat  engine.  Energy  is 
added  at  the  melting  point  and  the  maximum  rise  in  free  energy  dF  depends  on 
how  far  the  system  can  be  cooled  while  retaining  the  latent  heat  AH.  In  fact 
because  AF=<iH(dT/TJJ)  the  analogy  with  the  Carnot  efficiency  (AT/TM)  can  be 
made.  Like  the  heat  engine,  the  usable  stored  free  energy  is  reduced  by 
irreversible  processes  such  as  diffusion,  heat  flow,  and  interface  kinetics. 
Rapid  solidification  is  a  means  for  forstalling  solidification  and  achieving  a 
high  free  energy,  but  its  irreversible  aspects  dissipate  much  of  the  available 
free  energy.  Slow  cooling  with  controlled  nucleation  is  a  very  good  way  of 
producing  glass  or  highly  mecastable  crystalline  phases.  In  multicomponent 
solidification  slow  cooling  can  through  segregation,  produce  a  system  with 
much  free  energy  stored  in  the  wide  composition  range  and  in  the  eutectic 
phases  that  would  not  be  part  of  the  equilibrium  solid,  while  rapid  segregation¬ 
less  cooling  could  in  one  stage  produce  an  equilibrium  homogeneous  solid. 

The  controlled  ability  to  get  either  closer  to  equilibrium  or  further  from 
equilibrium  by  rapid  solidification  through  the  manipulation  of  the  kinetic 
and  thermodynamic  variables  must  surely  be  one  of  the  most  important  attributes 
of  this  new  processing  field. 

Multicomponent  Phase  Equilibria,  Partitionless  Phase  Changes 
and  Metastable  Hierarchies 

The  complication's  introduced  by  several  components  are  not  obvious  at  first 
glance.  Metastability  of  each  phase  is  again  common  and  limited  only  by 
spinodals  and  nucleation  and  growth  kinetics.  Thermodynamic  properties  of 
metastable  phases  are  well  behaved  and  often  smooch  extensions  of  behavior  in 
stable  ranges.  Equilibria  are  still  dictated  by  equality  of  chemical  potentials 
of  each  component  but  this  no  longer  implies  equal  free  energies.  A  phase 
diagram  having  as  axes  T,  P,  or  more  symmetrically  T,  uj,  v*2 . . . . iu. 

would  be  a  (N+l)-dimensional  analogue  of  the  single  component  diagram  where  N 
are  the  number  of  components.  Single  phases  would  occupy  (N+l)-dimensional 
hyper-volumes,  bounded  by  N-Jimcnsional  hypersurf aces  representing  two-phase 
equilibria,  which  intersect  at  (N-l)-dimensional  hypersurfaces  representing 
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three-phase  equilibria,  etc.  until  there  are  points  representing  (N+2)-phase 
equilibria.  Such  a  phase  diagram  is  useful  for  an  open  system  in  which  material 
reservoirs  keep  the  chemical  potentials  fixed  by  letting  compositions  change. 

In  a  few  situations  we  do  impose  open  system  conditions  for  some  of  the  com¬ 
ponents  by  e.g.,  fixing  the  partial  pressures  of  a  gas  containing  these  components, 
but  more  generally,  we  fix  the  composition  of  our  system,  which  leads  to 
profound  complications  in  the  thermodynamics.  The  open-system  metastable 
diagram  would  be  constructed  by  extension  of  all  two-phase  hypersurfaces  and 
would  divide  the  space  into  simplexes  in  each  of  which  a  hierachy  of  phase 
equilibria  would  exist. 

In  open  systems  the  entire  phase  diagram  is  occupied  by  single  phase 
hypervolumes  with  multiphase  equilibria  occupying  hypersurfaces  of  lower 
dimension  and  thus  of  zero  measure.  Unless  we  take  pains  to  fix  chemical 
potentials  and  temperature  to  lie  precisely  on  one  of  these  multiphase  surfaces, 
the  system  will  almost  always  be  single  phase.  If  it  is  multiphase  an 
infinitesimal  shift  in  conditions  off  the  surface  will  return  it  to  be  single 
phase.  The  open  system  metastable  phase  diagram  is  a  multidimensional  version 
of  a  one-component  phase  diagram  and  the  same  rules  including  the  hierarchial 
one  apply.  Indeed  Figures  2-4  could  serve  as  examples  of  a  two-dimensional 
section  of  a  multicomponent  phase  diagram. 

Major  changes  occur  when  we  shift  to  closed  systems  in  which  the  overall 
composition  rather  than  the  u's  are  fixed.  This  is  the  usual  situation  in 
rapid  solidification.  Phase  diagrams  for  closed  systems  have  compositions  as 
axes  and  contain  (N+l)-dimensional  hypervolumes  of  one-,  two-,  and  N-phase 
equilibria  with  (N+l)  and  (fH-2) -phase  equilibria  occupying  hypersurfaces  of 
zero  measure.  The  same  phases  remain  in  equilibrium  throughout  the  volume 
as  temperature,  pressure  and  composition  shift. 

A  common  graphic  method  of  depicting  the  multiphase  equilibria  is  the 
common  tangent  hyperplane  to  the  free  energy  hypersurfaces  for  the  various 
phases  at  a  given  temperature  and  pressure.  It  is  fully  equivalent  to  the 
condition  of  equal  chemical  potential  of  all  the  components,  and  the  range  of 
the  line  segment  between  tangent  points  readily  demonstrates  how  the  same  two 
phases  can  remain  in  equilbrium  over  a  wide  composition  range  (Fig.  5). 

The  stable  multiphase  equilibria  are  those  where  the  hyperplane  tangent  to 
the  free  energy  surface  of  the  equilibrium  phase  does  not  cross  the  free 
energy  surface  of  any  phase.  If  it  does  cross,  the  equilibrium  is  metastable 
with  respect  to  the  formation  of  the  phase. 

It  is  immediately  apparent  from  Figure  5  that  there  is  a  hierarchy  of 
metastable  phases  denoted  by  free  energy  curves  (or  hypersurfaces)  and  phase 
equilibria  denoted  by  tangent  line  segments  (or  hyperplanes).  The  boundaries 
to  a  given  hierarchy  are  given  first  of  all  by  the  tangent  points  on  the  free 
energy  curves.  These  tangents  map  out  the  stable  and  metastable  phase  diagram 
(Fig.  6). 

Three  other  features  mark  changes  in  the  hierarchy.  The  intersections  of 
free  energy  curves  (or  hypersurfaces)  define  the  T^  curve  (or  hypersurface) 
which  can  appear  on  phase  diagrams.  It  defines  a  restricted  equilibrium  for 
partitionless  transformation.  The  Tg  curves  do  not  denote  a  stable  or 
metastable  equilibrium.  They  do  pass  through  congruent  points  on  the  phase 
diagram  and  lie  within  the  corresponding  two-phase  field.  They  also  intersect 
at  Tq  triple  points  (or  hypercurves)  where  three  phases  have  the  same  free 
energy. 

T^  curves  in  Figure  6  mark  the  intersections  of  tangent  planes  and  therefore 
form  ruled  hypersurfaces  on  phase  diagrams.  T^  curves  mark  the  intersection 
of  a  surface  of  one  phase  with  a  plane  tangent  to  two  others. 
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Figure  5.  Two-component  free  energy  curves  for  various  phases  at 
constant  temperature.  The  left  and  right  portions  of  the  curves 
for  a  and  6  respectively,  and  the  tangent  between  them  represent 
the  stable  equilibrium.  Several  metastable  equilibria  are  depicted 
by  other  tangents.  The  crossing  of  free  energy  curves  represent 
points  on  the  T.  curves,  while  the  crossing  of  tangents  represent 
T.  curves.  Within  each  region  of  composition  separated  by  the 
vertical  lines  a  given  phase-equilibrium  hierarchy  holds. 


Figure  6.  A  binary-eutectic  diagram  at  constant  pressure  with 
metastable  extension  of  all  curves  and  the  T^  and  Tj  curves. 
Because  metastable  liquid  at  point  X  will  dissolve  sc.all  amounts 
of  the  stable  phase  S,  no  stability  hierarchy  based  on  phases 
per  se  can  exist. 


The  features  of  multicomponent  phase  diagrams  can  be  readily  imagined  to 
originate  from  free  energy  surfaces  of  various  phases  that  evolve  and  shift 
relative  to  each  other  with  changing  temperature  and  pressure.  It  is  seen 
that  the  metastable  continuation  of  a  stable  equilibrium  involves  only  the 
free  energy  curves  of  the  participating  phases  and  should  be  unaffected  by  the 
new  phase  whose  free  energy  surface  happened  to  move  through  the  tangent 
plane.  Thus  all  stable  equilibria  can  again  be  extended  to  depict  the 
metastable  equilibria  between  the  same  phases.  The  qualitative  evolution  of 
the  metastable  portion  of  Figure  6  can  be  understood  by  the  more  rapid  rise  in 
Figure  5  of  the  liquid  free  energy  relative  to  that  of  the  two  solids  with 
decreasing  temperature.  These  free  energy  curves  could  be  obtained  quantitatively 
by  a  number  of  methods  and  then  extrapolated  to  allow  determination  of  the 
metastable  diagram.  For  some  of  the  curves  it  seems  far  simpler  and  probably 
more  accurate  to  extrapolate  the  portions  of  a  measured  stable  (or  metastable) 
phase  diagram. 

As  for  the  single  component,  wherever  a  liquid  undergoes  a  glass  transition 
the  curves  will  be  continuous  witn  only  a  rapid  change  of  slope. 

It  should  be  immediately  pointed  out  that  no  simple  hierarchy  of  phases 
exists.  A  system  at  a  point  in  region  X  in  Figure  6  is  in  equilibrium  when  a 
and  6  are  present  as  the  stable  phases,  and  the  liquid  is  metastable.  Yet 
when  stable  8  is  brought  into  contact  with  liquid  in  X  under  conditions  where 
a  does  not  nucleate,  stable  S  will  dissolve  in  metastable  liquid.  Such  a 
reaction  in  which  a  stable  phase  dissolves  in  a  metastable  one  is  impossible 
in  single  component  systems,  but  quite  common  in  multicomponent  ones.  Only 
when  a  finally  nucleates  can  B  reappear.  This  disappearance  and  reappearance 
of  a  phase  in  spontaneous  processes  indicates  that  there  can  be  no  hierarchical 
listing  of  phases. 

This  can  also  be  indicated  by  listing  the  metastable  hierarchies  in  region 
X.  Because  X  lies  outside  the  6  +  1,  field,  no  equilibrium  between  6  and  L  is 
possible.  The  remaining  five  equilibria  are  ranked  as  follows:  single-phase 
B  is  highest  in  free  energy,  followed  by  single-phase  a.  single-phase  liquid, 
two-phase  a  +  L,  and  equilibrium  is  reached  with  two-phase  a  +  S.  In  region  W 
the  sequence  8  +  L  is  inserted  between  L  and  a  +  L.  In  such  metastable 
sequences  phases  can  appear,  disappear,  reappear  with  a  different  composition 
and  disappear  again.  Once  a  particular  phase  equilibrium  has  disappeared  it 
can  not  reappear  unless  the  system  is  reprocessed. 

The  inability  of  B  to  form  from  liquid  in  X  unless  a  has  also  appeared  is 
an  example  of  a  thermodynamically  required  precursor  reaction  (7) .  Such 
requirements  do  not  occur  in  single  component  systems,  because  any  of  the  more 
stable  phases  can  appear  at  any  time.  Consider  a  multicomponent  system  in 
metastable  equilibrium  of  one  or  more  (m)  phases.  The  composition  must  lie 
within  that  m-phase  field  on  the  metastable  phase  diagram.  The  addition  of  a 
new  phase  can  occur  only  if  the  composition  lies  also  beyond  where  the 
metastable  (m+l)-phase  field  first  intercepts  the  m-phase  field.  In  Figure  6, 

B  can  form  from  liquid  in  W  but  not  in  X,  while  a  can  form  in  either.  This 
places  a  strong  and  easily  formulated  constraint  on  the  phase-sequences  that 
can  occur  in  a  given  hierarchy  which  can  be  locally  applied  even  if  the  system 
has  only  reached  local  equilibrium.  - 

During  such  isothermal  sequences  the  free  energy  is  monotonically  decreasing, 
but  chemical  potentials  can  increase  as  well  as  decrease.  In  the  graphical 
construction  of  Figure  S,  the  chemical  potentials  of  the  two  components  are 
the  intercept  with  the  component  axes  of  a  tangent  to  the  single  phase  curve 
at  the  system  composition  or  of  the  common  tangent  in  a  multiphase  equilibrium. 
These  ups  and  downs  of  chemical  potentials  of  a  component  go  hand  in  hand  with 
the  appearance  or  disappearance  of  a  phase  enriched  in  that  component. 
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In  solidification,  thermodynamics  put  bounds  on  the  solid  composition  that 
can  form  from  liquid  given  the  composition  and  temperature  of  the  liquid  at 
the  interface  (1).  The  T.  curves  pertaining  to  liquid  with  a  given  solid  form 
an  upper  bound  to  the  solid  composition  that  can  form  from  liquid  of  any 
composition  at  that  temperature.  Because  the  T.  lies  between  the  liquidus  and 
solidus,  systems  in  which  the  extrapolated  solidus  appears  to  span  the  entire 
composition  range  are  good  candidates  for  segregationless  solidification  to 
that  particular  phase.  Systems  with  retrograde  solldi  (stable  or  metastable), 
even  where  the  corresponding  liquid!  sweep  across  the  composition  axis,  can  be 
shown  to  have  Tq  curves  which  are  bounded  in  composition.  Compositions  beyond 
this  bound,  which  is  approximately  the  liquidus  composition  at  the  retrograde 
temperature  (8),  can  not  undergo  segregationless  solidification  to  that 
crystalline  phase  at  any  temperature.  Because  crystallization  involving 
segregation  is  much  slower,  retrograde  systems  might  be  good  candidates  for 
glass  formers. 

Two  T^  curves  in  Figure  6  intersect  at  the  eutectic  point  Z.  The  eutectic 
horizontal  is  a  stable  T^,  the  other  moving  downward  in  temperature  is 
metastable  and  defines  bounds  to  conditions  where  it  is  impossible  for  one 
metastable  two-phase  equilibrium  to  follow  another.  A  T-  curve  (not  shown) 
also  goes  through  point  Z.  It  indicates  the  temperature  where  eutectic 
solidification  in  the  absence  of  a  proeutectic  reaction  first  becomes  thermo¬ 
dynamically  possible  although  unlikely. 

Because  of  the  segregation  accompanying  most  solidification  and  the  slowness 
of  diffusion  in  the  solid  state,  most  solidif ication  processes  yield  solids  in 
which  the  individual  phases  are  not  homogeneous.  The  hierarchy  diagrams  rank 
equilibria,  and  homogeneity  of  a  phase  is,  apart  from  minor  gravitational  and 
defect  segregation  phenomena,  a  necessary  condition  for  equilibrium  whether 
stable  or  mctastable.  These  segregation  effects  can  bring  the  system  into  a 
state  there  it  contains  phases  that  it  would  not  contain  in  any  metastable 
sequence. 

This  can  be  illustrated  by  an  example.  A  liquid  cooled  to  a  point  in 
region  Y  could  at  no  temperature  during  its  cooling  path  reach  any  equilibria 
ocher  than  a,  L,  L  +  a.  Yet  during  conventional  cooling,  segregation  will 
lead  to  the  formation  of  B  and  the  system  as  a  whole  will  not  reach  equilibrium 
until  all  the  6  has  disappeared.  The  hierarchy  diagrams  can  only  be  used  for 
conventional  solidification  on  a  local  equilibrium  basis.  They  are  more 
useful  for  the  solid  state  sequences  following  a  segregationless  solidification 
or  a  homogenized  solid  formed  by  rapid  solidification. 

The  usefulness  of  metastable  equilibrium  diagrams  lies  in. the  fact  that 
like  stable  diagrams  there  are  rules  for  their  construction  which  guide 
measurement  and  permit  our  experience  to  be  organized.  The  ability  to  extra¬ 
polate  and  interpolate  is  one  aspect  which  lets  us  make  rapid  strides  in 
sketching  in  the  main  feature  of  such  a  diagram.  When  metastable  equilibria 
is  reached,  thermodynamics  can  make  its  full  range  of  predictions.  For 
single  components  a  sequence  of  equilibria  are  likely.  The  possibility  of 
segregation  makes  it  difficult  for  us  to  bring  in  strong  predictions  about 
multicomponent  systems,  except  in  rapid  solidification  and  on  ^  local  level  in 
conventional  solidification. 

With  new  techniques  of  producing  phases  increasingly  far  from  stable 
equilibrium,  the  need  for  metastable  phase  diagrams  becomes  pressing.  While 
extrapolation  of  curves  from  stable  diagrams  has  been  long  used  for  estimation 
and  Tq  curves  have  been  in  use  for  diffusionless  transformation,  a  systematic 
approach  in  a  few  important  systems  should  prove  the  usefulness  of  the  ideas 
presented  here  in  charting  processes  far  from  stable  equilibrium. 
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RAPID  SOLIDIFICATION 


J.  W.  Cahn,  S.  R.  Coriell  and  W.  J.  Boettinger 
Rational  Bureau  of  Standards 
Washington,  D.C. 

Rapid  solidification  phenomena  are  described  in  terms  of  a 
hierarchi/  of  increasing  deviation  from  equilibrium.  Results 
of  moi-phological  stability  theoi'y  applied  to  silicon  regrowth 
indicates  that  factors  outside  of  conventional  constitutional 
supercooling  can  captain  the  observed  absence  of  lateral 
segregation.  A  model  for  interface  response  functions 
applicable  to  solute  trapping  is  presented  and  a  thermodynamic 
limit  to  the  amount  of  solute  trapping,  which  seems  to  be 
operating  in  some  systems,  is  proposed. 

X.  INTRODUCTION 

All  solidification  processes,  slow  or  rapid,  involve  some 
nonequilibrium  effects.  With  increasing  solidification  rates 
deviations  from  equilibrium  tend  to  become  more  extreme.  The 
extent  of  disequilibrium,  however,  is  nonuniform  and  different 
nonequi librium  phenomena  become  important  at  widely  differing 
rates.  We  shall  adopt  a  scheme  of  classifying  solidification 
effects  by  their  nonequilibrium  aspects. 

A,  Full  Diffusioral  Equilibrium  (Global  Equilibrium) 

At  the  very  slbwest  solidification  rates,  diffusional 
equilibration  can  occur.  In  this  state  of  equilibrium, 
materials  are  characterized  at  any  time  by  constant  temper¬ 
ature  and  chemical  potential  throughout  the  dimensions  L  of  a 
sample;  i.e.,  phases  are  homogeneous,  but  there  will  be  dis¬ 
continuities  in  composition  at  the  interfaces.  This  equilib¬ 
rium  requires  solid  state  diffusion  over  distances  of  the 
specimen  and  hence  geological  times  (L?/D.  where  D  is  the 
solute  diffusion  coefficient  in  the  sol id5phase).  sThe  ex¬ 
istence  and  compos i tions  of  homogeneous  phases  are  given  by 
stable  (or  mctastoble)  phase  diagrams.  This  kind  of  global 
equilibrium  will  not  be  approached  in  solidification  pro¬ 
cesses,  except  in  extremely  small  systems  such  as  ultrafine 
dispersions  of  droplets.  Surprisingly  because  the  time  scale 
for  nucleation  can  also  be  of  geological  times,  metastable 
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phases  can  occur  in  this  regime  (!)•  In  classical  thermo¬ 
dynamics  metastable  equilibrium  is  a  full  equilibrium  and 
indistinguishable  from  stable  equilibrium  except  that  one  or 
more  stable  phases  are  prevented  from  nucleation  (2). 

B.  Local  Equilibrium 

This  next  range  of  phenomena  is  characterized  by  the 
absence  of  diffusional  equilibration  in  either  solid  or  melt. 

Each  element  of  the  system  is  assumed  to  have  reached  a  state 
where  it  can  be  adequately  described  by  the  local  state 
variables  of  temperature,  pressure  (stress),  and  composition. 

The  interface  is  also  in  local  equilibrium.  Temperature  and 
chemical  potentials  of  each  species  can  be  assumed  continuous 
across  the  interface,  and  thus  the  limiting  compositions  of 
each  phase  at  the  interface  will  be  given  by  a  knowledge  of 
the  interface  temperature  and  the  equilibrium  phase  diagram, 
with  appropriate  corrections  for  the  effect  of  interface 
curvature.  These  limiting  compositions  then  form  boundary 
conditions  to  the  transport  equations  within  each  phase  and  a 
wide  range  of  solidification  processes  can  be  modelled  with 
these  assumptions  (3),  Within  each  phase  local  equilibrium 
can  be  established  in  times  of  the  order  of  the  diffusional 
jump  frequency.  At  the  interface  depending  on  orientation  and 
defect  structure  there  are  a  number  of  relaxation  processes 
with  widely  differing  time  scales  (e.g.,  nucleation  of  a  new 
layer  of  crystal),  that  can  lead  to  early  breakdown  of  local 
equilibrium  there  (4). 

C.  Deviations  from  Local  Equilibriwn 

Effects  due  to  deviations  from  local  equilibrium  at  the 
interface  already  become  obvious  in  some  systems  at  very  small 
undercoolings  and  very  small  interface  velocities.  For  ex¬ 
ample,  dislocation- free  crystals  in  contact  with  slightly 
superheated  or  supercooled  melts  show  no  measurable  growth 
rates.  With  increasing  interface  disequilibrium,  growth  or 
melting  will  begin  (4,5).  While  the  facetted  form  can  be  an 
equilibrium  form  in  contact  with  a  fluid  of  uniform  temper¬ 
ature  and  chemical  potential,  a  growing  facetted  crystal 
rarely  is  at  lccal  equilibrium,  since  melt  isotherms  or  satura¬ 
tion  contours  are  not  likely  to  follow  a  facetted  shape  (4). 

At  equilibrium,  segregation  of  solutes  between  solution  and 
crystal  is  independent  of  interface  orientation.  Thus  the 
observation  of  compositional  discontinuities  in  Czochralski 
grown  crystal  between  regions  that  grew  behind  facets  (core) 
and  regions  that  grew  behind  smoothly  curved  (off-cere)  liquid 
melt  interfaces  indicates  that  nonequilibrium  interface 
effects  are  frequently  observed  at  conventional  crystallization 
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rates  (6).  The  prediction  of  the  extent  of  nonequi librium 
interface  phenomena  has  thus  been  an  important  goal  which  has 
received  renewed  attention  with  the  investigation  of  rapid 
solidification. 

This  goal  may  be  summarized  by  response  functions  (7) 
which  for  example  give  local  growth  velocities  and  solid 
compositions  at  the  interface  as  a  function  of  conditions  at 
the  interface.  These  response  functions  replace  local 
equilibrium  conditions  at  the  interface  and  become  the  boundary 
conditions  to  the  transport  equations  in  each  phase.  ' 

While  it  is  clear  that  many  rapid  solidification  effects 
can  only  be  understood  if  the  concept  of  local  equilibrium  is 
discarded,  some  predictions  based  on  local  equilibrium  have 
not  been  fully  explored.  In  the  present  paper  we  describe  a) 
results  of- morphological  stability  theory  relevant  to  silicon 
regrowth  at  high  velocities  using  local  equilibrium  and  b) 
results  of  a  model  for  the  interface  response  functions 
applicable  to  solute  trapping. 

II.  MICROSEGREGATION  AND  STABILITY  AT  THE  SOLIDIFICATION 

INTERFACE  * 

During  the  freezing  of  alloys,  a  planar  solid-liquid 
interface  can  become  unstable.  The  solute  rejected  by  the 
growing  crystal  can  lead  to  a  thin  layer  of  undercooled-1 iquid 
surrounding  the  crystal  which  then  makes  the  shape  unstable 
with  respect  to  dendrites  that  could  protrude  into  this  thin 
layer.  This  leads  to  a  cellular  or  dendritic  interfaces  and 
causes  solute  microsegregation.  The  classical  theory  of 
constitutional  supercooling  (3,8)  predicts  increasing  in¬ 
stability  with  increasing  velocity  (for  constant  liquid 
temperature  gradient)  and  would  predict  cellular  or  dendritic 
solidification  for  most  rapid  solidification  experiments. 

There  are,  however,  several  assumptions  that  may  not  be  valid 
for  rapid  solidification.  (1)  Local  equilibrium  at  the  inter¬ 
face  is  assumed  to  describe  quantitatively  the  amount  of 
solute  rejected  by  the  growing  crystal.  (2)  Capillarity  and 
latent  heat  evolution  which  could  stabilize  the  interface  are 
neglected.  (3)  The  theory  indicates  a  layer  of  liquid  cooler 
.  than  its  local  melting  point,  but  does  not  show  how  the 
growing  crystal  breaks  through  to  this  layer. 

Beginning  with  the  work  of  Mullins  and  Sekerka  (9-12),  a 
method  of  analysing  morphological  instabilities  has  evolved 
which  properly  solves  the  transport  equations  assuming  local 
equilibrium  at  the  interface  or  making  some  allowance  for 
deviation  from  equilibrium.  The  principal  result  is  that 
there  is  far  greater  stability  during  rapid  solidification 
than  predicted  by  consti lutional  supercooling  even  if  local 
equilibrium  is  maintained.  The  transition  between  planar  and 


nonplanar  growth  during  the  directional  solidification  of 
binary  alloys  can  be  calculated  from  linear  morphological 
stability  theory  assuming  local  equilibrium.  For  constant 
velocity  V  solidification  in  the  z-direction,  the  perturbed 
solid-liquid  interface  is  of  the  form 

z  a  5  exp  (ot  +  iuxx  +  iuyO, 

where  6  exp  (ot)  is  the  amplitude  and  w  and  u  are  spatial 
frequencies  of  the  perturbed  interface.  If  the  real  part  of  a 
is  positive  for  any  perturbation  (for  any  values  of  u  and 
u  ),  the  interface  amplitude  increases  exponentially  Oith  time 
tyand  the  planar  interface  is  unstable.  If  the  real  part  of  a 
is  negative  for  all  possible  perturbations,  the  interface  is 
stable.  The. value  of  o  can  be  determined  by  solving  the  heat 
flow  and  diffusion  equations  with  appropriate  boundary  con¬ 
ditions.  Assuming  local  equilibrium  at  an  isotropic  solid- 
liquid  interface,  Sekerka  (13)  found 

o  =  V{-kLGL(aL-V/KL)  r  ksGs(os+V/<s)  -  2kTMru,2a  + 

2TotG  a(a-V/D)/(a-pV/D)}/{L  V  +  2kkla/(a-pV/D)}  (1) 

c  V  c 

with 

a  *  (V/2D)  +  £(V/2D}2  +  w2  +  o/D]**, 

*  (V/2<^)  +  [(V/2k^)2  +  a»2  +  o/k^]*5 
os  =  -(V/2k$)  +  C(Y/2<S)2  +  u2  +  o/Ksj* 
o  «  (ks°s  +  kLaL)/(2k), 

where  <,  and  <  are  thermal  diffusivities  of  liquid  and  solid, 
respectively,  k,  and  k  are  thermal  conductivities  of  liquid 
and  solid,  respectively,  0  is  the  liquid  diffusion  constant 
(diffusion  in  the  solid  is  neglected),  u  *  (w  2  +  u  2)’*,  k  = 
(ks+k,)/2,  G,  and  G,.  are  the  unperturbed  .temperatur^  gradients 
in  the  liquid  and  solid,  respectively,  G  is  the  unperturbed 
solute  gradient  i.e.,  G  *  Vc^(k-1 ) / 0 k , kc i s  the  equilibrium 
distribution  coefficient,  coi  Ts  the  bulk  solute  concentration, 
Ty  is  the  melting  point  of  the  planar  interface  in  the  absence 
of  solute,  r  the  ratio  of  the  solid-liquid  surface  energy,  > , 
and  the  latent  heat  of  fusion  per  unit  volume,  L  ,  m  is  the 
slope  of  the  liquidus  line  on  the  phase  diagram, vand  p  =  1-k. 

For  a  given  alloy  and  specified  values  of  the  processing 
conditions  we  can  calculate  numerically  u  as  a  function  of 
'’(>••).  from  the  above  equation.  Alternatively,  wo  set  o  •-  0 
and  calculate  the  bulk  concentration  c  as  a. function  of 
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c— (<*•) •  The  minimum  value  of  c^w)  gives  the  stability- 
instability  demarcation. 

The  sign  of  o  is  determined  by  the  sign  of  the  numerator 
in  the  above  equation  since  the  denominator  is  always  positive. 
The  numerator  consists  of  four  terms  proportional  to  G.  ,  G  , 
r,  and  G  ,  respectively.  For  the  usual  case  of  positive 
temperature  gradients,  the  first  two  terms  are  negative  re¬ 
flecting  the  stabilizing  influence  of  the  temperature  field. 

The  surface  tension  term  is  always  stabilizing.  Instability 
occurs  when  the  destabilizing  solute  gradient  term  is  suffi¬ 
ciently  large  to  overcome  the  stabilizing  influence  of  the 
temperature  gradients  and  surface  tension. 

It  is  usually  an  excellent  approximation  to  neglect  the 
thermal  diffusivity  terms  in  the  above  equations,  i.e.,  let 
ie.->  •»  and  -*■  ®.  In  this  limit,  Sekerka  (14)  has  shown  that 
the  stability  criterion  is 

2kLGl  +  VLv  >  <ks+kL)mGcS(A,k)‘  (2) 

where  S(A,k)  is  a  function  of  the  dimensionless  variables  k 
and  A=(kTMrV2)/(D2mG  ).  The  function  S(A,k)  can  be  calculated 
by  solving  a  cubic  equation;  among  its  properties  are  that  0  <_ 
S(A,k)<.  1,  S(0,k)  =  1  and  S(A  >  l,k)  =  0.  The  stability- 
instability  demarcation  for  silicon  contair.g  arsenic  and 
antimony  for  G.  =  2G0  K/cm  is  indicated  by  the  curve  in  Figure 
1.  L 

The  following  values  were  used  for  the  calculations; 

kL  *  70  Jin'1  s'1  iC1,  ks=  22.0m'1 s'1  k"1 , 

D  *  3xl0_8m2/s,  k  =  0.3  (As). .k  »  0.02  (Sb), 

Ly  *  4.56  GJ/m3,  T^r  *  1.3x1  O'7  Km, 

*  3xl0’5m2/s,  «s=  9.4xl0“6m2/s, 

m  =  -440K/at.  fract.(As),  m  =  -370  K/at.  fract.  (Sb). 

For  a  given  velocity  V  for  solute  concentrations  below  the 
curve  the  planar  solid-liquid  interface  is  stable.  The 
straight  line  with  negative  slope  in  the  figure  is  the  con¬ 
stitutional  supercooling  criterion,  mG  =  G.  or  c  «  G,  0k/ 
[V(k-l)m],  for  interface  stability.  This  fol lows“*frosh  the 
above  equation  in  the  limit  S(A,k)  -►  1,  k  -  k.  ,  and  2k.  G,  >> 

VI  .  For  Si,  it  appears  that  even  at  modest  velocities,  the 
preceding  inequality,  i.e.,  G.  >>3.26xl0’;V(sK/m' ) ,  is  not  a 
good  approximation,  and  the  latent  heat  term  in  the  stability 
criterion  can  be  important'.  The  straight  line  with  positive 
slope  is  the  Mull ins-Sckerka  absolute  stability  criterion  (9), 
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FIGURE  1.  The 
critical  concentration 
of  a)  >ts  and  b)  Sb 
above  which  interface 
instability  occurs  as 
a  function  of  the 
interface  velocity  V 
of  directional  solidi¬ 
fication  of  silicon  for 
cr  temperature  gradient 
in  the  liquid  of  2C0 
K/ctn.  The .  curve  is 
based  on  morphological 
stability  theory  while 
the  lines  with  negative 
and  positive  slope 
correspond  to  the  con¬ 
stitutional  super¬ 
cooling  and  absolute 
stability  criteria 3 
respectively . 
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As1  or  mG  *  kT..r(V/D)2  or  =  k2T..rV/[D{k-l )m].  For  small 
velocities,  the  constitutional  supercooling  criterion  is  a 
good  approximation  to  the  morphological  stability  calculations. 
For  very  large  velocities  the  absolute  stability  criterion  is 
a  good  approximation.  Mullins  and  Sekcrka  have  shown  that  the 
exact  curve  will  always  lie  above  the  absolute  stability  line. 
Thus*  for  solute  concentrations  below  those  given  by  the 
absolute  stability  equation,  the  interface  is  stable.  Whereas 
the  constitutional  supercooling  criterion  depends  linearly  on 
the  temperature  gradient  in  the  liquid,  G,  ,  the  absolute 
stability  criterion  is  independent  of  G.  \'although  we  require 
that  k  G  +  k.  G.  >  0).  In  Figure  2  we  show  the  wavelength, 

\  -  2it/w,  at  the  onset  of  instability  as  a  function  of  velocity. 
The  wavelength  decreases  monotonically  with  increasing  velocity. 
Although  we  show  results  for  V  =  10  m/s,  the  theory  is  no 
longer  valid  at  such  high  velocities  since  the  diffusion 
boundary  layer  length,  D/V,  is  only  a  few  atomic  distances. 

The  theory  has  been  extended  to  include  nonlinear  effects 
(12),  some  limited  interface  anisotropy  and  deviations  from 
local  equilibrium  (15,16). 


FIGURE  2.  The  wavelength. 

\  at  the  onset  of  instability 
during  directional  solidifi¬ 
cation  of  silicon  containing 
arsenic  or  antinony  as  a 
f motion  of  interface  velocity 
V. 
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III.  NONEQUILIBRIUM  AT  THE  INTERFACE  • 

Crystal  growth  or  dissolution  can  only  occur  with  some 
deviation  from  equilibrium  at  the  interface.  The  response  of 
the  interface  to  a  given  deviation  from  equilibrium  has  been 
called  the  response  function  (7).  Theories  attempt  to  predict 
for  single  component  systems  a  relationship  between  interface 
undercooling,  orientation,  crystal  perfection,  components  of 
interface  curvature  and  the  local  interface  velocity.  For  a 
binary  system  an  additional  relation  is  needed  to  link  the 
solid  and  liquid  compositions  to  the  above  variables. 

For  the  single  component  systems  theories  have  concen¬ 
trated  on  mechanisms  of  crystal  growth  which  broadly  fall  into 
two  classes:  continuous  additions  to  the  crystal  or  lateral 
spreading  of  new  crystal  layers.  Crystal  perfection  only 
matters  for  the  latter  class  of  mechanisms,  in  which  case 
imperfect  crystals  (containing  certain  types  of  dislocations 
or  twins)  could  grow  much  closer  to  equilibrium  than  perfect 
crystals.  A  rule  based  on  entropy  of  fusion  and  packing  on 
planes  parallel  to  the  interface  has  considerable  success  in 
approximately  classifying  those  systems  in  which  lateral 
growth  is  likely  (17).  It  has  been  predicted  that  at  high 
undercoolings  many  of  these  systems  change  mechanisms  and  grow 
by  the  continuous  mechanism  (4). 

For  binary  systems  an  early  kinetic  theory  (18)  predicted 
for  the  continuous  growth  mechanism  the  growth  rate  and 
partitioning  of' solute  in  terms  of  four  rate  constants 
describing  the  frequencies  with  which  the  two  kinds  of  atoms 
jump  from  melt  to  solid  and  v.ice  versa.  An  examination  of  its 
predictions  showed  that  this  theory  predicted  that  the  solid 
composition  could  never  exceed  the  equilibrium  composition  of 
solid  in  contact  with  liquid  at  the  interface  temperature 
(19),  To  exceed  the  solid  solubility  maximum  at  a  eutectic, 
solidification  would  have  to  occur  at  temperatures  where  the 
extrapolated  metastable  solidus  curve  exceeded  the  solid 
composition.  A  test  of  this  theory  required  either  knowing 
the  interface  temperature  during  solidification  or  choosing  a 
system  with  a  maximum  in  the  solidus,  a  system  with  retrograde 
solid  solubility.  Several  such  systems  have  now  been  examined 
and  the  finding  of  solid  solutions  in  excess  of  the  maximum 
equilibrium  solidus  conclusively  shows  that  this  kind  of 
theory  is  inadequate  (19-21). 

Such  a  kinetic  theory  leads  to  the  result  that  each  com¬ 
ponent  must  experience  a  decrease  in  chemical  potential  at  the 
interface,  if  that  component  is  to  jump  into  the  solid.  The 
experience  with  the  retrograde  systems  proves  that  the  minor 
component  can  increase  its  chemical  potential  upon  solidifica¬ 
tion,  a  condition  termed  solute  trapping  (7,19-21).  This 
condition  does  not  violate  thermodynamics  and  is  easily 
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modelled  if  one  assumes  that  the  solute  is  passively  carried 
into  the  solid  by  the  solidification,  but  can  respond  to  this 
indignity  by  diffusing  back  into  the  liquid  (22).  Baker 
worked  out  several  such  models  including  a  fairly  interesting 
one  that  examined  a  wide  variety  of  parameters  and  showed  how 
sensitive  the  predictions  of  such  models  are  to  minor  adjust¬ 
ments. 

This  model  is  shown  in  Figure  3.  The  interface  is  a 
moving  snowplow  for  the  solute,  whose  energy  is  higher  in  the 
solid.  Following  a  suggestion  by  Chernov  (23)  that  trapping 
might  occur  by  having  the  interface  be  a  state  of  low  solute 
energy  to  which  the  solute  voluntarily  rushes  only  to  be 
(passively)  buried  by  the  addition  of  the  next  layer,  Baker 
permitted  a  variable  solute  energy  at  the  interface.  He  also 
allowed  for  diffusion  in  both  phases.  The  resulting  con¬ 
tinuous  steady  state  diffusion  equation  with  drift  term  when  E 
is  a  function  of  z  is  readily  integrated  (24,25)  to  give  the 
composition  profiles  as  a  functional  of  velocity,  E(z)  and 
0(z).  Taking  the  ratio  of  the  limiting  compositions  of  liquid 


FIGURE  3.  Baker's  model  for  solute  redistribution  allowed 
for  an  arbitrarily  varying  diffusion  coefficient  and  potential 
energy  for  solute  near  the  moving  interface.  The  particular 
form  depicted  above  leads  to  the  redistribution  given  in 
equation  (3)  and  depicted  in  Figure  4. 


and  solid  at  the  interface  leads  to  the  following  expression 

(22). 

(CL/Cs)int  *  {exp(-a-e)  +  expt-fOCkg/kj-expt-ci)]/ 


[1  +  ln(ke/ki)/a]  +ke[l-exp{-[5)/ki]/[l+lnki/B]}/ke  (3) 

where  . 

ke  =  exp[(EL-Es)/kT] 
ki  -  exp  [{E^E^/kT] 

are  the  equilibrium  segregation  coefficients  between  liquid 
and  solid  or  interface  respectively  and 

a  =  «5V/DS 

e  =  sv/dl 

are  reciprocals  of  diffusion  distances  (in  units  of  the  bound¬ 
ary  thickness  5)  in  solid  and  liquid  respectively.  Figure  4 
shows  the  result  for  k  =0.1  and  D  /D^  =  10”J. 


FIGURE  4.  Partitioning  at  interface  as  a  function  of 
imposed  growth  rate  for  the  model  depicted  in  Figure  5  (from 
(22). 


Let  us  first  examine  curves  3  and  4  for  which  E.  is  E  and 
E.  respectively.  For  these  special  conditions  the  energys 
curve  in  Figure  3  has  a  single  ramp  which  occurs  in  a  region 
of  liquid  or  solid  diffusion  respectively.  Comparison  of  the 
two  shows  the  effect  of  shifts  of  the  energy  ramp  relative  to 
the  discontinuity  in  diffusivity.  When  the  solute  is  in  the 
energy  ramp  it  experiences  a  drift  force.  In  curve  4  this 
drift  is  exerted  when  the  diffusivity  is  much  lower  than  it  is 
for  the  model  of  curve  3.  Trapping  occurs  at  lower  velocities. 
For  curves  1  and  2  the  solute  is  first  repelled  from  the 
interface.  Segregation  over  a  wide  range  of  velocities  is 
stronger  than  equilibrium  until  at  the  highest  velocities 
trapping  finally  occurs.  Chernov's  ideas  are  modelled  in 
curves  5  and  6.  The  solute  is  attracted  to  the  interface  and 
then  buried'  in  the  solid.  Trapping  occurs  at  quite  low 
velocities,  but  in  the  intermediate  velocity  range  segregation 
to  the  solid  is  so  strong  that  the  effect  overshoots.  Baker 
concluded  that  Chernov's  ideas  were  not  necessary  to  explain 
trapping.  He  also  concluded  that  very  accurate  modelling  of 
interface  structure,  kinetics  and  energetics  would  be  neces¬ 
sary  to  obtain  valid  predictions  about  solute  trapping  except 
at  the  very  highest  velocities.  Curves  3  and  4  are  quite 
similar  to  the  model  presented  by  Jackson,  ct  al.  (26)  at  this 
conference. 

The  thermodynamic  questions  raised  by  solute  trapping  were 
explored  extensively  by  Baker  and  Cahn  (7).  Of  particular  in¬ 
terest  is  the  Tq  curve  which  is  the  intersection  of  the  solid 
and  liquid  freeuenergy  plotted  on  a  phase  diagram.  It  not 
only  marks  a  bound  for  diffusionless  solidification  in  which 
solid  and  liquid  can  have  the  same  composition;  it  also  marks 
the  bound  on  solid  compositions  which  can  form  from  liquid  of 
any  composition  at  that  temperature.  Where  thermodynamic  data 
are  accurate  enough  to  predict  metastable  extensions  of  free 
energy  functions,  TQ  curves  can  be  obtained.  For  the  retro¬ 
grade  systems  a  simple  estimate  can  be  constructed  for  the 
maximum  composition  that  the  TQ  curve  reaches  at  0°K.  No 
solid  of  this  phase  in  excess  of  this  maximum  can  form  from 
the  liquid  no  matter  how  rapidly  cooled,  unless  local  equilib¬ 
rium  ceases  for  the  liquid  phase  itself. 

This  estimate  is  obtained  from  a  rigorous  relationship 
(27,2S)  for  the  slope  of  the  solidus  curve. 


dT  _  (ycs)T(^Gs/jcs») _ 

<lcs  Hl(cl)-C(1-cl)hu(cs)+clH2s(cs)] 


(4) 


whore  H.  (c,  )_is  the  molar  enthalpy  of  liquid  of  composition  c, 
and  ll.  and  !!.,s  are  the  partial  molar  enthalpies  of  the  two  u 
components  in’" the  solid  at  composition  c  and  G  is  the  molar 
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Gibbs  free  energy  of  the  solid.  At  a  given  temperature  the 
denominator  has  a  simple  geometric  interpretation  on  an  enthalpy 
vs.  composition  plot  for  the  liquid  and  solid  phases  (wnic  -  is 
relatively  independent  of  temperature).  The  denominator  is 
the  vertical  distance  at  ttie  liquidus  composition  for  that 
temperature  between  the  liquid  enthalpy  curve  and  the  tangent 
drawn  to  the  solid  enthalpy  curve  at  the  solidus  composition. 
Hence,  at  the  retrograde  temperature,  where  the  slope  of  the 
solidus  is  infinite,  the  denominator  is  zero  and  the  tangent 
crosses  the  liquid  enthalpy  curve  at  the' liquidus  composition 
(q)  in  Figure  5.  At  low  temperatures  where  enthalpy  is  the 
dominant  term  in  the  free  energy,  the  TQ  composition  can  be 
approximated  by  the  intersection  of  theenthalpy  curves.  In 
the  absence  of  any  information  about  the  curvature  of  the 
solid  enthalpy  to  determine  the  intersection  (p)  of  the  en¬ 
thalpy  curves,  an  approximate  limit  on  the  T^. composition  can 
be  obtained  by  the  intersection  of  the  previously  mentioned 
tangent  and  the  liquid  enthalpy  curve  (q).  This  intersection 
is  the  liquidus  composition  at  the  retrograde  temperature. 

The  T0  curve  is  an  absolute  limit  to  solute  trapping  and  its 
exacvlocation  is  important  for  many  rapid  solidification 
predictions  (31).  Such  a  limit  seems  to  have  been  observed  by 


—  —  Rctrogrodi  T«mp*iotur« 


Ape*0«»mc»«  Bound* 
to  Tg  Curve 


Composition 


FIGURE  5.  Us  ins  the 
rigorous  relatzonsmp  sc. 

(4)  betveen  solidus  slope 
and  enthalpy  curves,  the 
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the  intersection  (q).  This 
composition  approximates  pz. 
(p)  the  intersection  of  the 
liquid  aid  solid  enzhaipy 
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of  the  Tq  carpositicn  at 
absolute  zero. 
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Whi  to  et  al .  (20).  In  Table  I  we  give  the  maximum  equilibrium 
solid  solubility  at  the  retrograde  c  ,  the  liquid  composition 
at  the  retrograde  c,  ’  ,  or  where  it  isn't  known  an  estimate 
c  r/kn,  and  the  maximum  solid  composition  that  has  been 
observed  to  date.  Except  for  Si (As),  c  cX  is  less  than 
c. 1 (or  c  r/kfl),  consistent  with  the  estimated  thermodynamic 
bound  given  above.  Every  one  of  these  systems  shows  solute 
trapping  but  T0  is  a  limit  to  how -far  this  can  be  pushed 
regardless  of  the  solidification  rate.  These  systems  are 
approaching  our  estimate  of  this  limit.  It  will  be  worthwhile 
to  increase  solidification  rates  still  further  to  test  the 
existence  of  this  limit.  When  calorimetric  or  free  energy 
data  for  these  systems  become  available  a  more  precise 
estimate  of  the  composition  at  the  crossing  of  the  enthalpy 
curves,  which  is  the  rigorous  basis  of  the  thermodynamic 
limit,  should  be  substituted  for  c^r. 
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Table  I.  Comparison  of  the  Liquidits  Composition  at  the 
Retrograde  Tarperature  with  the  Maximum  Solid 
Solubility  Found  in  Rapid  Solidification 


Solvent 

(solute) 

c  r 
s 

c  r/k 
cs  /K0 

^  max 
c 

<* 

Observed 

Comments 

Zn(Cd) 

0.025 

.22(29) 

.10(19) 

mole 

f motions 

Si  (As) 

l.SxlO21 

SxlO21 

6x1 021 
(20) 

atoms/acj  k. 
from  (30j 

Si(Sb) 

7x1 01 9 

3x1 022 

1.3xl021 

(20) 

atoms /oa , kn 
from  (30) 

Si(Ca) 

IS 

4.5x10 

6xl021 

A  C  1*20 

9.  SxlO 
(20) 

atoms/cc, kQ 
fivm  (20) 

Si (In) 

5x10 17 

2x10 21 

"0 

lxlQ-U 

(20) 

atoms /ec. k  . 
from  (30)  0 

Si (Si) 

t-xic27 

1x1 021 

nO 

•1x10 

(20) 

mono  tactic 
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